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Abstract 

A recent development in finite temperature field theory, the so-called 
Braaten-Pisarski method, and its application to properties of a quark-gluon 
plasma, possibly formed in relativistic heavy ion collisions, are reviewed. 
In particular parton damping rates, the energy loss of energetic partons, 
thermalization times, viscosity, and production and damping rates of hard 
photons are discussed. 
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Chapter 1 
Introduction 



There are two different methods based on the fundamental theory of strong interac- 
tions, i.e. QCD, for a theoretical description of the quark-gluon plasma (QGP) and 
predictions of possible signatures, namely lattice and perturbative QCD. The advan- 
tage of lattice QCD (Creutz, 1983) is the non-perturbative nature of this method, 
rendering the entire regime of values of the coupling constants and temperatures 
accessible. However, lattice simulations are numerically very elaborate. Thus only 
relatively small lattice sizes with typically 16 3 X 8 sites (Gottlieb et al., 1993) can be 
treated up to now. Hence the extrapolation of the results to the continuum case is not 
established completely. Moreover, it is not possible to describe dynamical variables on 
the lattice, such as signatures of the QGP based on particle production. Furthermore, 
neither a finite quark chemical potential nor pre-equilibrium effects can be considered 
on a lattice. On the other hand, the equation of state and the phase transition can 
be studied within lattice calculations (Petersson, 1991). A strong variation in the 
energy density, corresponding to a hadron gas below the critical temperature and to 
a QGP above, respectively, has been found. The critical temperature resulting from 
these investigations is given by T c = 150 - 200 MeV corresponding to a critical energy 
density of e c = 1 - 3 GeV/fm 3 . Above about twice to four times the critical temper- 
ature the lattice equation of state approaches the one of an ideal parton gas. The 
determination of the order of the phase transition is more involved. Recent results 
taking into account quarks dynamically indicate rather a smooth transition (cross 
over) than a first or second order phase transition (Brown et al., 1990). 

Perturbative QCD requires a small value of the strong coupling constant and 
works only above the critical point. Estimates of the effective temperature dependent 
coupling constant between T c and 2T C lead to a s = g 2 /47r = 0.2 - 0.5. (The lower 
value can be motivated by renormalization scheme arguments at finite temperature 
(Huang and Lissia, 1994).) With increasing temperature the coupling constant should 
decrease - presumably logarithmically - according to asymptotic freedom. Perturba- 
tive QCD can be extended to finite temperature (and chemical potential) (Kapusta, 
1989) for investigating properties of a thermalized QGP, e.g. dynamical signatures. 
Small deviations from the equilibrium, e.g. thermalization times and transport coeffi- 
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cients, can be considered using linear response theory (Heinz, Kajantie, and Toimela, 
1987) or the Kubo formalism (Zubarev, 1974). 

However, applying perturbation theory to gauge theories at finite temperature 
leads to serious problems. Many physical quantities, computed perturbatively, suffer 
from infrared singularities and gauge dependence. For example, different results for 
the damping rate of a gluon at zero momentum were found in different gauges, even 
negative ones in covariant gauges (Lopez, Parikh, and Siemens, 1985), correspond- 
ing to an unphysical instability of the QGP. This so-called plasmon puzzle was the 
starting point for the development of an improved, effective perturbation theory by 
Braaten and Pisarski (1990a). Pisarski (1989b) discovered that results obtained in 
naive perturbation theories can be incomplete in the order of the coupling constant 
since higher order diagrams can contribute to lower order in the expansion in the cou- 
pling constant. Braaten and Pisarski (1990a) succeeded, by distinguishing between 
hard momenta of the order of T and soft momenta of the order of gT, to isolate the 
relevant diagrams for QCD. These so-called "Hard Thermal Loops" (HTL) have to 
be resummed in a geometrical series (see also Frenkel and Taylor (1990).) In this 
way effective propagators and vertices are constructed, which can be used as in or- 
dinary perturbation theory leading to results that are complete to leading order in 
the coupling constant. This effective perturbation theory yields gauge independent 
results for physical quantities (Braaten and Pisarski, 1990c) and in the case of the 
gluon damping rate a positive value (Braaten and Pisarski, 1990b). At the same time 
medium effects are taken into account due to the resummation, e.g. Debye screening 
caused by the color charges of the QGP, removing or at least weakening infrared 
divergences. Quantities that are logarithmically infrared divergent to leading order, 
using bare Green's functions, are finite applying the Braaten-Pisarski method. Many 
interesting quantities of the QGP belong to this class, such as the energy loss of en- 
ergetic partons in the QGP (Braaten and Thoma, 1991b), the viscosity of the QGP 
(Thoma, 1991a), and the production rate of hard photons in the QGP (Kapusta, 
Lichard, and Seibert, 1991; Baier et al., 1992). Quantities that are quadratically in- 
frared divergent in naive perturbation theory, e.g. damping rates of hard partons, are 
logarithmically infrared divergent using effective Green's functions (Pisarski, 1989b; 
Thoma and Gyulassy, 1991). The reason for this behavior is due to the fact that the 
effective gluon propagator contains screening of static chromoelectric fields (Debye 
screening) but no static chromomagnetic screening. A consistent consideration of the 
latter screening probably requires non-perturbative methods, such as a self consistent 
resummation based on a mean field approximation. So far there are only preliminary 
attempts for developing such methods for QCD (Kalashnikov, 1992; Thoma, 1992b). 

Another problem of the Braaten-Pisarski method relies on the assumption of 
the weak coupling limit, 3 < 1, for distinguishing between hard and soft momenta. 
To what extent the results obtained under this assumption can be extrapolated to 
realistic values of the coupling constant, g = 1.5 - 2.5, has to be studied in each 
single case. Anyway, the Braaten-Pisarski method means an enormous improvement 
compared to the naive perturbation theory at finite temperature and takes account 
of important medium effects of the QGP. 
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Last but not least, investigations of field theoretic methods for thermal systems 
are of fundamental interest. For this problem requires the combination of the three 
basic methods of modern physics, i.e. statistical physics, theory of relativity, and 
quantum mechanics. The application of relativistic quantum statistics or thermal field 
theory, however, is restricted to relativistic systems, such as plasmas of relativistic 
particles, which appear in nature only under extreme astrophysical conditions, e.g. 
supernovae explosions, neutron stars, and the Big Bang. Owing to the possibility of 
creating a QGP in relativistic heavy ion collisions the development of field theoretic 
methods at finite temperature has been intensified within the last few years greatly. 

The aim of the present report is the introduction to the Braaten-Pisarski method 
and to review its application to interesting problems of the QGP in ultrarelativistic 
heavy ion collisions. In order to provide a useful introduction to this topic for non- 
experts who want to apply this method, we will try to keep the discussion on an 
informal level as simple as possible. (An alternative, more formal introduction to 
thermal field theories and the Braaten-Pisarski method has been presented by Altherr 
(1993).) Some basic ideas of thermal field theory and perturbation theory at finite 
temperature, relevant for the following, will be discussed in chapter 2. The Braaten- 
Pisarski method will be introduced in chapter 3, where we will avoid complicated and 
formal derivations as far as possible. For this purpose the basic ideas are exemplified 
using the scalar field theory. The remaining part of this report is devoted to the use 
of the Braaten-Pisarski method for interesting properties of the QGP, demonstrating 
at the same time the technical details of the method and its application. In chapter 
4 the damping rates of particles in the QGP are treated, which are the starting point 
for the energy loss of partons in the QGP (chapter 5) and the transport rates (chapter 
6), from which thermalization times and the viscosity of the QGP follow. In chapter 
7 the damping and production of hard photons will be considered. The conclusions 
including the limitations of the Braaten-Pisarski method and possible extensions are 
given in chapter 8. In appendix A the notations used throughout this report are 
summarized. In appendix B the photon self energy in the HTL approximation will 
be calculated explicitly. 
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Chapter 2 

Perturbation theory at finite 
temperature 

Our goal is the computation of cross sections, life times, etc. of relativistic particles 
at high temperatures. High temperature means a temperature T that is significantly 
higher than the masses m of the thermal particles, T 3> m, leading to a relativistic 
heat bath due to particle production, such as a QGP consisting of it-, J-quarks and 
gluons. For the description of such a system we will adopt perturbative methods 
similar to the zero temperature case, i.e. Feynman diagrams and rules for computing 
amplitudes. 

2.1 Imaginary time formalism 

Green's functions are the central objects of perturbation theory. The most important 
Green's function is the propagator. How does it look like at finite temperature? For 
simplicity sake, we first consider the scalar field theory 1 . At zero temperature the 
bare propagator is given by the vacuum expectation value of the time ordered product 
of two fields at different space-time points (Mandl and Shaw, f 984) 

iA T = o (x-y) = (0\T{<P(x)<P(y)}\0) , (2.1) 

where |0) denotes the ground state of the non-interacting theory. At finite tempera- 
ture the expectation value in (2.1) has to be replaced by a thermal one, e.g. using a 
canonical partition function Z = tr[exp(—/3H)] with f3 = 1/T and the Hamiltonian 
H: 

i A(x - y) = \ 5>TO(*My)}l»> c-^» , (2.2) 

where E n and \n) are the eigenvalues and eigenstates of the non-interacting Hamilto- 
nian. Expressing the field by creation and annihilation operators in the four dimen- 
sional momentum space, K = (k 0} k), the expectation value in (2.2) can be evaluated, 

1 The free Lagrangian reads Co = (<9 |[1 </><9' I </> — m 2 <f> 2 )/2. 
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for instance for x > y 0} 

iA(x-y) = j 7^7^ {[! + n B M] e~ U<{x ~ y) + n B (to k ) e iK ^} (2.3) 
with k = LOk = Vk 2 + m 2 and the Bose- Einstein distribution 

»*M = • ( 2 - 4 ) 

Eq.(2.3) describes the propagation of a scalar particle in a heat bath from x to y } 
where besides the spontaneous emission in the vacuum at x there is also induced 
emission and absorption described by the terms containing nfl(o;t). 

As in the zero temperature case (Mandl and Shaw, 1984), we want to express 
the propagator by a four dimensional integral over K. A possible representation can 
be obtained in the so-called imaginary time formalism, which we will use here, since 
the Braaten-Pisarski method has been derived within this formalism. Owing to the 
periodicity of the thermal propagator in imaginary time r = it, A(r) = A(r + f3) 
(Kapusta, 1989), the Fourier integral for k is replaced by a Fourier series 

£ liSfio^^' (2 - 5) 

n= — oo V / 

where the four momentum K contains now discrete energies k = 2TrinT with integers 
n (Matsubara, 1955). The equivalence of (2.5) and (2.3) follows from the relation 
(Kapusta, 1989) 

1 f °° j,. If(ko) + f(-ko)] 



T J2 f(ko = 2TrinT) = — / dk { - 

n=-oo 27TZ J-joo 



1 pioo-\-e 

+ 77- / . dk [f(k ) + f(-k )] n B (k ) , (2.6) 

Z7TZ J — ioo-\-e 

which holds for arbitrary functions f(k ) having no pole on the imaginary axis. 

Compared to the zero temperature Feynman rules, the propagator in momentum 
representation is now given by A(A^) = 1/(K 2 — m 2 ) with the Matsubara frequencies 
k = 2TrinT and the loop integration by iTJ2k I d 3 k/(2Tr) 3 . In the case of fermions 
the propagator is anti-periodic, implying k = (2n + l)iirT . 

As simplest example we consider the tadpole self energy of the </> 4 -theory 2 at 
finite temperature shown in fig.2.1. Using standard Feynman rules (see appendix A) 
modified by the finite temperature prescription given above, we find 

d 3 k 1 



The zero temperature part vanishes after renormalization. The finite temperature 
part is ultraviolet finite due to the distribution function and can be evaluated exactly 
in the limit of a vanishing mass yielding the simple result 

III = g 2 T 2 . (2.8) 



An interaction term —g <f> is added to the free Lagrangian. 
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Figure 2.1: Tadpole self energy of a scalar particle with </> 4 interaction. 

2.2 Saclay method 

Sums over k in loops, containing two or more propagators, can be performed most 
conveniently by the so-called Saclay-method (Pisarski, 1988). Here the propagator is 
presented by 

A(K) = - fdre k » r A(r,cu fc ) , (2.9) 
Jo 

where 

A(r,co fc ) = -T^e- fcoT A(/0 (2.10) 

= 7T^ {[1 + "sWl e"^ T + n B (to k ) e^ T \ . 

The last equation follows from using (2.6). In this way the summation over k in 
expressions containing A(if) is trivial, because k shows up only in the exponential 
function in (2.9). 

In the case of a fermion propagator the Saclay representation reads 

S(P) = (P^ + M)A(P) (2.11) 

with 



and 



A 



A(P) = - f dr e P0T A(r, E p ) (2.12) 
Jo 

(r,£ F ) = ^ {[1 - n F {E r )] e~ E ^ - n F (E p ) e E ^} , (2.13) 



where E p = \/p 2 + M 2 and 

ME P ) = (2.14) 

denotes the Fermi-Dirac distribution. It is straightforward to take into account a 
finite chemical potential (see e.g. (Vija and Thoma, 1994)) in addition. 
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Figure 2.2: Photon self energy to lowest order. 



2.3 Self energies 



The self energy of a gauge boson in the one-loop approximation is an important exam- 
ple for the application of the perturbation theory at finite temperature. In general it 
is not possible to derive an analytic formula due to distribution functions, appearing 
under the integrals. In the high temperature limit, however, i.e. for temperatures 
much larger than the momentum of the gauge boson, the integrations can be done ex- 
actly. The gluon self energy in the high temperature limit has been found by Klimov 
(1982) and independently by Weldon (1982a). 

For simplicity sake we will hrst treat the photon self energy (fig. 2. 2). We consider 
a thermal system of relativistic electrons, positrons and photons with a net charge 
zero, i.e. a vanishing chemical potential. We will follow the derivation by Braaten and 
Pisarski (1990a), which relies on the restriction to hard loop momenta k,k ~T. As 
discussed in appendix B this hard thermal loop (HTL) approximation is equivalent 
to the high temperature limit. 

The transversality of the polarisation tensor, P fi H fll/ (P) = 0, implies that only two 
components are independent, for which we choose the longitudinal and the transverse 
components 



The explicit computation using the methods developed in the last sections is shown 
in appendix B for the longitudinal part. The transverse part follows analogously. The 
leading term of the high temperature expansion is given by 



where m 7 = eT/3 can be regarded as a "thermal photon mass", generated by the 
interaction of the photons with the thermal electrons and positrons. 
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Figure 2.3: Gluon self energy to lowest order. 



The following remarks are of interest: 

1. The next term of the high temperature expansion can also be calculated exactly 
(Toimela, 1986). It is proportional to g 2 T times a function of p and p, which has the 
dimension of an energy. 

2. The photon self energy is momentum independent in contrast to the tadpole of 
the </> 4 -theory (2.8), which is also proportional to g 2 T 2 . 

3. For pi < p 2 } (2.16) has an imaginary part 



The resulting damping is called Landau damping (Pisarski, 1988) and is essential for 
quantities like the energy loss or the viscosity, as will be discussed later on. 

4. The self energy (2.16) reduces to Hl(po = 0) = —3m 2 and Ht(po = 0) = in 
the static limit p — *■ 0. As discussed in the next section, the static limit of the self 
energy is related to the screening of static fields, and m D = 3m 2 is the Debye mass, 
leading to the screening of static electric fields, whereas there is no screening of static 
magnetic fields within the HTL approximation (Kajantie and Kapusta, 1985). This 
result is not surprising since an electron-positron gas contains electric charges but no 
magnetic monopoles. 

5. The leading term of the high temperature expansion can be derived also from the 
Vlasov equation, i.e. a transport equation without collision term considering a mean 
held. The only non-classical quantity, entering thereby, is the Fermi distribution for 
the electrons and positrons. The photon self energy is obtained now via the dielectric 
constants (Elze and Heinz, 1989; Mrowczyhski, 1990; Thoma and Gyulassy, 1991) 
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Figure 2.4: Electron self energy to lowest order. 



which follow from the solution of the the Vlasov equation for the distribution func- 
tions assuming a small deviation from equilibrium (Lifshitz and Pitaevskii, 1981). 
The coincidence of the self energy, found in the perturbative high temperature ap- 
proximation, with the one of the Vlasov equation is caused by the equivalence of the 
high temperature limit, T — > oo, and the classical limit, % — > 0. As a matter of fact, 
already Silin (1960) found the result (2.16) in this way. 

Next we turn to the gluon self energy. For this purpose the diagrams of fig. 2. 3 
have to be considered. Because of the gluon propagator a gauge has to be chosen. 
Surprisingly, the leading term of the high temperature expansion differs from (2.16) 
just by replacing the thermal photon mass m 7 by a thermal gluon mass (Klimov, 



where Nj indicates the number of thermal quark flavors. In contrast to the photon 
self energy (fig. 2. 2), which is of course gauge independent, gauge independence holds 
for the gluon self energy only for the leading term of the high temperature expansion 
(Heinz, Kajantie, and Toimela, 1987). The coincidence of the QED and QCD polari- 
sation tensor can be understood from the Vlasov equation, where the only difference 
between QED and QCD originates from the distribution functions (Elze and Heinz, 



The electron and quark self energies in the HTL limit can be derived analogously 
by starting from the diagram in fig. 2. 4 (Klimov 1982, Weldon 1982b). Neglecting the 
bare fermion mass M it reads 



1982; Weldon, 1982a) 





1989). 





with 



1 

ip 2 



[tr(P^E)-p tr( 7o E)] , 





where the traces over the 7 matrices are given by 



HP^ E) 



4 ml 
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Figure 2.5: Dyson-Schwinger equation defining the effective propagator of a scalar 
particle. 

tr( 7o £) = 2ml-\n^^- (2.22) 
P Po ~ P 

with the thermal fermion masses m\ = e 2 T 2 /8 in the QED and m\ = g 2 T 2 /6 in the 
QCD case. 

2.4 Effective propagators and dispersion relations 

First we consider the ^-theory with a vanishing mass. The effective propagator is 
defined by the Dyson-Schwinger equation in fig. 2. 5, where the tadpole (2.8) shows up 
in the r.h.s. The effective propagator is obtained by iteration leading to a geometrical 
series, in which the tadpole is resummed. The result of this resummation is given by 

A-(A) = K2 1 2 , (2.23) 

where m 2 th = TIi = g 2 T 2 is the thermal mass. The dispersion relation following from 
the pole of the effective propagator, 

u(k) = \Jk 2 + m 2 h , (2.24) 

describes the propagation of a collective mode with mass m^, caused by the interac- 
tion with the heat bath. 

Analogously to the ^-theory we construct an effective photon propagator. It is 
convenient to choose the Coulomb gauge, d{ A{ = 0, for this purpose. 

The bare propagator in Coulomb gauge is given by 



A»(A) = -, 



f 

Am (A) = 0, 

D 13 (K) = 1- U 3 - k ~§) • (2-25) 
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Figure 2.6: Dyson-Schwinger equation defining the effective photon propagator. 



Here the Doo-term describes the propagation of a longitudinal photon mediating the 
electric interaction (Coulomb held). The transverse part By corresponds to the prop- 
agation of a transverse photon and the mediation of the magnetic interaction, respec- 
tively. 

The effective photon propagator is defined by the Dyson-Schwinger equation 
shown in fig. 2. 6 yielding the longitudinal and transverse contributions 

1 



DUK) 



D* T (K 



k 2 



1 



k 2 - n T (K) 



(2.26) 



where we assumed the same Lorentz structure as for the bare propagator, i.e. D^ = 
D% D* oi = and D% = D* T (% - k t k 3 /k 2 ). ' 

Eq.(2.26) represents the most general form of the exact photon propagator in 
Coulomb gauge (Heinz, Kajantie, and Toimela, 1987). Now we assume the high 
temperature expressions (2.16) for the self energy II ^ and II x in (2.26). The dispersion 
relation of the collective photon mode in the relativistic electron-positron plasma 
follows from the poles of the effective propagator (Klimov, 1982) 



P + Zm* 1 



id id + k 

— In 

2k to 



k 







to 



k 2 



6 2 U 

2 m -F 



1 



1 - -r -r In 



id + k 



to 



k 



. (2.27) 



Before we are going to solve these equations, we will consider the static limit, 
to — > 0, of the effective photon propagator. The longitudinal part reduces to D* L (u — > 



0) = l/(k + 3m^), the transverse to D\ 



[to 



0) 



-Ijk . Thus the long range 



electric interaction is shielded in the infrared limit by the Debye mass m 2 D = 3m 2 , 
originated by the presence of electric charges in the medium, as opposed to the mag- 
netic interaction. 

The equations (2.27) cannot be solved analytically. In the limits, k — > and 
k — > oo, however, there are simple solutions (Pisarski, 1989a): 



oo L {k -> 0) 



3 k 2 

m-y H 

7 10 m, 



13 



co/m 




Figure 2.7: Dispersion relations of collective photons and gluons. 



U>i,(k — > Oo) 

uj T (k -> 0) 



k, 

3 P 
m 1 + - — , 



5 m~ 



oo 



A; . 



(2.28) 



The numerical solution of (2.27) is depicted in fig. 2. 7. Let us discuss these results: 
f. In the limit k — > oo the dispersion relation of a bare photon, u = k } is recovered. In 
the limit A; — > the dispersion relation describes a non-relativistic particle with mass 
Comparing with a non-relativistic plasma (Jackson, 1975) one observes that this 



m 7 . 

particle corresponds to a collective mode called plasmon in the the longitudinal case. 
The thermal photon mass m 7 is identical to the plasma frequency. 
2. The dispersion relation lies above the light cone, u > k } for all values of k. 
Therefore the imaginary part of the self energy does not contribute and the plasma 
oscillations are not damped in the high temperature approximation; i.e., there is no 
Landau damping, causing the dissipation of energy from the plasma wave into the 
heat bath (Lifshitz and Pitaevskii, 1981). 

Analogously to the collective photon modes in a QED plasma there are collective 
gluon modes in a QGP, which simply result by replacing m 7 in (2.27) and (2.28) by 
/// .. of (2.19). 

Besides collective bosonic modes also collective fermionic modes exist in a rel- 
ativistic plasma following from the poles of the effective fermion propagator. It is 
convenient to adopt the helicity representation for the effective propagator of an elec- 
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tron or quark (Braaten, Pisarski, and Yuan, 1990): 



where 

D±(P) = - Po ±p + -^ [±tr(P^, E) - {± Po - p) tr( 7o £)] . (2.30) 

The traces containing the fermion self energy are given by (2.22). The dispersion 
relations are given by the zeros of D±(P) (Braaten, Pisarski, and Yuan, 1990). It is 
interesting to note that the mode corresponding to D_(P), called plasmino (Braaten, 
1992), has no zero temperature counterpart, as it is also the case for the plasmon. 



2.5 Gluon damping rate 

The gluon damping rate caused a controversy in the middle of the eighties about the 
consistency of perturbative QCD at finite temperature (plasmon puzzle), triggering 
investigations that led to a consistent formulation of perturbation theory within the 
Braaten-Pisarski resummation technique (Pisarski, 1989b). Therefore we will sketch 
the status of the gluon damping before the invention of the Braaten-Pisarski method. 

The gluon damping rate is defined as the imaginary part of the dispersion relation, 
*y(k) = —Imu(k). Here we restrict ourselves to the damping of collective modes 
at zero momentum in a pure gluon gas. Assuming no overdamping, 7 <C Reuj } the 
damping rates of the transverse and longitudinal modes are given by (Heinz, Kajantie, 
and Toimela, 1987; Braaten and Pisarski, 1990b) 



7x(0) 



1 



2m, 



Im Ht(lo = m g , k = 0) 



7L(0) = -^lim 



Im II l (uj = m g , k) 



(2.311 



Since there is no direction specified in the plasma, the dielectric tensor, 



1 i \ kikj 



(2.32) 



reduces to e(cu) Sij implying e£,(u;,& = 0) = ex(cu, k = 0) (Lifshitz and Pitaevskii, 
1981). Then (2.18) together with (2.31) yields 7t (0) = 7 L (0). 

As we have seen, there is no gluon damping associated with the leading term of 
the high temperature expansion; i.e., the dispersion relation has no imaginary part, 
because u > k is true for all k and the imaginary part of the gluon self energy in the 
high temperature limit vanishes for to 2 > k 2 . The next to leading term of the self 
energy in the high temperature expansion, on the other hand, contains an imaginary 
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Figure 2.8: Damping of a collective gluon according to (2.33) and (2.34) caused by 
absorption of a thermal gluon and the decay in two thermal gluons, respectively. 

part, leading to damping. This imaginary part of the longitudinal contribution in the 
limit, k «C tt>, which is of interest here, reads (Heinz, Kajantie, and Toimela, 1987) 

i p 

ImH L {k<.Lo) = g 2 T — . (2.33) 

47T U 

This result has been derived within the temporal and Coulomb gauge, respec- 
tively. As mentioned already above, the non-leading terms of the high temperature 
expansion are gauge dependent. The gluon damping rate at zero momentum can be 
written as 

7L,x(0) = a 9 -^- , (2.34) 

87T 

where a = 1 holds in temporal and Coulomb gauge. Using covariant gauges, however, 
negative results have been found, e.g. a = —5 in the case of the Feynman gauge 
(Lopez, Parikh, and Siemens, 1985). Negative damping rates would correspond to 
an instability of the QGP. Moreover, observables like damping rates should be gauge 
independent. These unphysical results led to speculations about the failure of per- 
turbative QCD at hnite temperature. However, Pisarski (1989b) noticed that the 
perturbation theory, as it was used here, is incomplete; i.e., diagrams contributing to 
the same order, g 2 } in (2.34) have been neglected. The modification of perturbation 
theory, for taking consistently into account all contributions to the same order, is the 
topic of the next chapter. 

Before we will come to this, the mechanism for the gluon damping according to 
(2.34) shall be discussed. The physical meaning of the imaginary part of the gluon self 
energy on the mass shell can be illuminated according to an extension of the Cutkosky 
cutting rules to hnite temperature (Weldon, 1983), by cutting the internal lines of 
the diagram. We consider the diagram corresponding to polarisation by gluon pair 
creation. Owing to the high temperature expansion the external momentum has to 
be small compared to the temperature, while the internal momentum is of the order 
of the temperature (see appendix B); i.e., we consider the polarisation of a collective 
mode, indicated by a blob in fig. 2. 8, caused by thermal gluons. Cutting the diagram 
we obtain to possible processes: (1) damping of a collective mode by absorption of 
a thermal gluon and (2) decay into two thermal gluons. The hrst process, requiring 
space like momenta, K < 0, is forbidden because of u > k for the collective gluon. 
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The second process is possible only if the energy of the final gluons is equal to the 
one of the collective gluon. Since the decay takes place in the plasma, the final gluons 
have then to be collective modes as well. As a matter of fact, this reflects already the 
inconsistency in the derivation of (2.34). The true damping mechanism to order g 2 T 
is based on higher processes, as it will be discussed in section 4.2. 
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Chapter 3 

The Braaten-Pisarski method 



As we have seen in the last chapter, studying the example of the gluon damping 
rate at zero momentum, the use of bare Green's functions at finite temperature may 
lead to unphysical results. This inconsistency of the naive perturbation theory relies 
on the fact that higher order diagrams, i.e. multi-loop diagrams, can contribute to 
a lower order in the coupling constant at finite temperature than naively expected. 
These diagrams are the HTL diagrams, of which we encountered already the polar- 
isation tensor, (2.16), and the fermion self energy, (2.20) to (2.22). They have to 
be resummed into effective Green's functions, which will be used then just like bare 
Green's functions in naive perturbation theory. This effective perturbation theory, 
developed by Braaten and Pisarski (1990a) leads to consistent results that are com- 
plete in the order of the coupling constant. In the following this idea will be discussed 
in detail. 

3.1 </) 4 -theory 

In order to exemplify the program, sketched above, we will start with the </> 4 -theory. 
Neglecting the mass of the scalar particle - we are interested in gauge theories with 
massless bosons - we found 111 = g 2 T 2 , (2.8), for the tadpole self energy. Now 
we want to go beyond the one-loop self energy. In the two-loop approximation the 
diagram shown in fig. 3.1 arises, which is proportional to g 4 and infrared divergent 
in the massless case, because there are two propagators in the lower loop. Naive 
expectation suggests that the exact self energy can be written as II = g 2 T 2 [1 + 0(g 2 )] } 
providing a correction of the order g 2 to the tadpole. 

In order to check this conjecture, we consider the tadpole containing the effective 
propagator (2.23) instead of the bare one (fig. 3. 2). This effective tadpole II* consists 
of the sum of all daisy diagrams. Apart from the bare tadpole 111 all these diagrams 
are infrared divergent. In spite of containing an infinite sum of diagrams the effective 
tadpole can be computed easily by simply replacing the mass m in (2.7) by the thermal 
mass m t h = gT . 
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Figure 3.1: Contribution to the two-loop sell energy ol a scalar particle. 



First we will consider the ultraviolet hnite contribution. From (2.7) we get 



n%n = -g 2 / dk — n B (n k ) 

1 7T JO ilk 



(3.1) 



with 0| = k 2 -\-m 2 h . This integral cannot be evaluated exactly, but it can be expanded 
into a series lor small g. For this purpose we decompose the integral by introducing a 
separation scale /.• ' . restricted by gT <C k* <C T, and make use ol fe, m t h <C T in the 
hrst integral and ol m t h <C k in the second. In this way one folds (Kapusta, 1989) 



This result also holds it the ultraviolet divergent part is added, since it is ol order g 
alter renormalization. 

The result (3.2) shows some interesting properties: 

1. The effective tadpole II* is infrared hnite, although it contains (intmitely many) 
infrared divergent diagrams. 

2. The correction to the bare tadpole IT is ol the order ol g instead ol g 2 } as naively 
expected. 

3. The correction is not a multiple ol the coupling constant g 2 ol the theory and thus 
not a perturbative result. 

4. An expansion in the number ol loops is not equivalent to one in the coupling 
constant, as diagrams with an arbitrary number ol loops contribute to g 3 . 

The reason lor this surprising behavior ol the sell energy is the use ol the effective 
propagator, containing an infrared regulator mth = gT in the denominator, caused by 
the resummation ol intmitely many diagrams (see hg.2.5), implying a hnite result lor 
IP. Furthermore, the order in the coupling constant is reduced, since II* is sensitive to 
small momenta and energies; i.e., small momenta and energies contribute significantly 
to the integral in (3.1) in such a way that the coupling g in the denominator ol the 
integrand in (3.1) becomes important. This behavior, exemplified in the case ol <f> 4 - 
theory, is typical lor held theories at hnite temperature, as long as the bare masses 
ol the particles are zero or negligible compared to the temperature. 
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Figure 3.2: Scalar tadpole containing the effective propagator. 



These observations motivate the following recipe, consisting of three steps, for 
massless held theories at hnite temperature: 

1. First, diagrams proportional to g 2 T 2 have to be isolated, in our example the 
tadpole 111 in (2.8). 

2. Next, effective Green's functions have to be constructed by resummation of these 
diagrams, here the effective propagator A* of (2.23). 

3. Finally, quantities such as self energies - here the effective tadpole II* (3.2) - have 
to be calculated, using an effective perturbation theory, in which, depending on the 
situation, bare or effective Green's functions are used as in ordinary perturbation 
theory. If only momenta and energies of the Green's functions are required which 
are hard, i.e. of the order of the temperature or larger, bare Green's functions are 
sufficient, as can be seen for example from the effective scalar propagator A*, in 
which the effective mass m 2 h = g 2 T 2 can be neg lected if K 2 ~ T 2 . (For an exact 
argumentation it is convenient to adopt the Euclidean metric, for which K 2 = k^ + k 2 
holds. Then K 2 is always hard, if k and (or) k are hard (Braaten and Pisarski, 
1990a).) If, however, all momenta and energies of the Green's functions, necessary 
for the calculation, are soft, i.e. proportional to gT } effective Green's functions have 
to be used. This is the case for the scalar self energy beyond the leading order g 2 T 2 } 
since small momenta are important in the loop integration of the effective tadpole. 

It is essential to emphasize that this method is based on the weak coupling as- 
sumption, </ --v., 1 . which allows for a distinction between hard and soft momenta. 

The proceeding, proposed here, is the idea of the effective perturbation theory 
developed by Braaten and Pisarski (1990a). In the next section it will be extended to 
gauge theories, which are of course more involved than the </> 4 -theory. Besides the <f> 4 - 
theory the Yukawa theory (Thoma, 1994c) and scalar electrodynamics (Kraemmer, 
Rebhan, and Schulz, 1994) have been studied as simple examples for applying the 
Braaten-Pisarski method. 

3.2 Gauge theories 

We follow mainly the arguments given by Braaten and Pisarski (1990a). The starting 
point is the fact that naive perturbation theory, i.e. the exclusive use of bare Green's 
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Figure 3.3: HTL correction to the bare gluon propagator. 



functions, does not reproduce all contributions to leading order in g for certain quan- 
tities, such as the gluon damping rate. Therefore one has to ask the question, which 
diagrams are as important as the tree diagrams defining the bare Green's functions. 
In other words, which diagrams lead to effective Green's functions that contribute to 
the same order in g as the bare ones? The answer to this question depends on the 
external momenta of these diagrams. If these momenta are hard, P ~ T, there are no 
such diagrams besides the bare Green's functions. If, however, the external momenta 
are soft, P ~ gT then diagrams exist from which effective Green's functions can be 
constructed that must not be neglected compared to the bare ones. These diagrams 
are the HTL. 

As an example we consider the gluon propagator of fig. 3. 3, where the gluon self 
energy in the HTL approximation proportional to g 2 T 2 has been inserted. If the 
external momentum is soft, this propagator, consisting of two bare propagators and 
the HTL polarisation tensor, is as important as the bare one. If, however, P is hard, 
it gives contributions reduced by a factor of g 2 . It should be noted that the internal 
momenta are hard, since only hard momenta contribute to the HTL gluon self energy 
(see appendix B). 

Now we want to proceed according to our recipe, proposed in the last section, by 
isolating the HTL diagrams first. In the case of the ^-theory the only HTL diagram 
is the momentum independent tadpole, Hi = g 2 T 2 . In gauge theories one encounters 
two complications. First, the HTL are momentum dependent, e.g. the polarisation 
tensor (2.16). Secondly, the gauge symmetry leads to relations between the Green's 
functions (Ward identities), causing effective vertices in addition to effective propa- 
gators. For instance, the relation between the gluon self energy and the three-gluon 
vertex is given by 



Hence there are HTL corrections to the three-gluon vertex, shown in fig. 3. 4, which 
are as important as the bare vertex. Analogously we expect HTL corrections to the 
four-gluon vertex and the quark-gluon vertex. Braaten and Pisarski (1990a) showed 
that HTL corrections arise for all n-gluon amplitudes and n-gluon-plus-one-quark- 
pair amplitudes. In Coulomb or Feynman gauge these corrections are represented by 
the one-loop diagrams with hard loop momenta, containing only three-gluon vertices 
and quark-gluon vertices (with the exception of the tadpole in fig. 2. 3), as shown in 
fig. 3. 5. (Since this method leads to gauge independent results (Braaten and Pisarski, 



(K - a% r^ p (A, a") = n^A) - w(k') . 



(3.3) 
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Figure 3.4: HTL corrections to the three-gluon vertex. 



1990c), it is possible to choose any gauge.) These diagrams have to be taken into 
account in the corresponding amplitudes only if all external legs are soft. This can be 
proven by power counting in g similarly to fig. 3. 3 without calculating the diagrams 
explicitly. 

The computation of these HTL is very tedious due to their complicated momen- 
tum and energy dependence. However, it can be simplified using Ward identities 
(Braaten and Pisarski, 1990d). Furthermore, it is possible to derive the HTL from 
a generating functional, which can be obtained by adding the following term to the 
QCD Lagrangian (Braaten and Pisarski, 1992a) 



where m g and m q denote the effective gluon and quark masses defined in (2.19) and 
(2.22), D the covariant derivative, and the angle brackets averaging over the directions 
of P p leading to the logarithmic factors in the HTL (see e.g. (2.16)). This effective 
Lagrangian shows a couple of interesting properties, which are the topic of current 
investigations (Brandt et a!., 1993). 

Moreover, the HTL can be derived from a classical kinetic theory similar to the 
polarisation tensor, where the only non-classical inputs are the Bose and Fermi dis- 
tributions (Blaizot and Iancu, 1993; Kelly et al., 1994). 

According to the next step of our recipe effective Green's functions will be con- 
structed from the HTL. In the case of the effective gluon or quark propagator this 
happens via resummation within a Dyson-Schwinger equation, which has be done 
already in section 2.4 (see (2.27), (2.29) and (2.30)). (Inserting the HTL self en- 
ergy only once as in fig. 3. 3 is not sufficient, since repeated insertion of Hhtl in the 
case of soft external momenta results obviously in effective propagators which are as 
important as the bare ones.) 




(3.4) 
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Figure 3.5: All HTL corrections of QCD in Coulomb- or Feynman gauge. 



For producing the effective vertices, the HTL corrections have to be added simply 
to the bare vertices. In fig. 3. 6 all effective Green's functions up to the four-point 
amplitude are listed. Here also an effective two-quark-two-gluon vertex appears which 
has no bare contribution. This vertex plays a role for the plasmino damping (see 
section 4.2). 

The last step of our recipe consists of the use of the effective Green's functions 
within an effective perturbation theory for computing certain quantities, e.g. prop- 
erties of the QGP. When do we have to use effective Green's functions and when are 
bare ones sufficient? As soon as soft momenta become important in the diagrams 
which contribute to the quantity under consideration, i.e., as soon as the quantity 
becomes sensitive to the momentum scale gT effective Green's functions have to be 
considered. However, only if all legs attached to the Green's functions are soft, it is 
necessary to use effective Green's functions. For example, the gluon self energy in 
the HTL approximation (2.16) requires only bare Green's functions since only hard 
momenta in the loop contribute. Hence the construction of effective Green's functions 
from bare ones is consistent. Otherwise, a self consistent resummation would have 
been necessary. 

For calculating the next to leading order of the gluon self energy proportional 
to g 3 } however, effective propagators and vertices have to be utilized (Rebhan, 1993; 
Schulz, 1994). In the case of the gluon damping rate they have to be taken into 
account already to leading order (see section 4.2). As a matter of fact, each quantity 
becomes sensitive to momenta of the order of gT from a certain order on and has 
then to be treated applying the Braaten-Pisarski method. From which order on this 
occurs, has to be studied in each single case. 

In this way one obtains a systematic expansion in the coupling constant and, as 
shown by Braaten and Pisarski (I990d), gauge independent results for amplitudes 
on the mass shell. (The latter statement has been doubted recently, since seem- 
ingly gauge dependent results have been found for the quark damping rate (Baier, 
Kunstatter, and Schiff, 1992). This problem will be discussed in the next chapter.) 

In conclusion, the effective perturbation theory of Braaten and Pisarski means a 
crucial improvement compared to the naive perturbation theory at finite temperature. 
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gure 3.6: Effective Green's functions of QCD up to the four-point amplitude. 



For physical quantities it provides gauge independent results that are complete in the 
order of the coupling constant. Furthermore screening effects due to the self energies 
in the denominator of the effective propagators, e.g. Debye screening, are taken 
into account, thus improving the infrared behavior significantly. Many interesting 
observables which are infrared divergent within the naive perturbation theory are 
finite using the Braaten-Pisarski method. 

Nevertheless, the effective perturbation theory is not the ultimate solution of all 
problems of thermal field theory, since it suffers from two serious disadvantages. First, 
it is based on the assumption j < 1, which is never fulfilled in the case of the QGP. 
Typical values for situation, met in ultrarelativistic heavy ion collisions, are given 
by g = 1.5 - 2.5, corresponding to a s = g 2 /4:ir = 0.2 - 0.5. If one assumes that 
the effective temperature dependent coupling constant decreases logarithmically with 
increasing temperature, g < 0.5 is obtained at a temperature of the order of the Planck 
mass (Miiller, 1993). To what extent the extrapolation of the results, obtained via the 
Braaten-Pisarski method, to realistic values of the coupling constant can be justified, 
is debatable and has to be investigated in each single case. On the other hand, the 
requirement j < 1 might be too restrictive (Braaten, 1995), because it depends on 
the definition of the coupling constant, which might contain numerical factors as 
for instance in the case of the ^-theory, where one often chooses </ 2 /4! instead of 
g 2 . A criterion for the validity of the effective perturbation theory, which does not 
depend on such a convention for defining the coupling constant, might be derived by 
comparing the inverse of the bare gluon propagator in the imaginary time formalism 
k 2 + {2'wnT) 2 with the static longitudinal propagator k 2 + m 2 D . Demanding that the 
soft momentum scale is significantly smaller than the hard one, m 2 D <C (27rT) 2 , implies 
a s <C 7r instead of <C 1. Indeed, the applications of the Braaten-Pisarski method 
to the quantities discussed in chapters 5-7 are valid as long as g ~ 1 holds. 

In any case, however, the Braaten-Pisarski method takes account of important 
properties of the QGP, e.g. medium effects, at least qualitatively. Moreover, it pro- 
vides an improvement of the naive perturbation theory which gives results for example 
for the equation of state above about two to four times the critical temperature in 
agreement with lattice calculations (Petersson, 1991). Anyway, the problem of the 
magnitude of the coupling constant is not a conceptual problem of the method but 
one of the application. 

The second problem, on the other hand, is a problem of principle, inherent in the 
way the resummation has been performed. The Braaten-Pisarski method considers 
no magnetic screening since the transverse part of the HTL polarisation tensor (2.16) 
vanishes in the static limit, p — > 0, as discussed already in section 2.3. Therefore 
it can happen that certain quantities are infrared divergent at leading order even 
when calculated with effective Green's functions. This occurs if the quantity becomes 
sensitive to momenta of the order of g 2 T } i.e., if momenta of the order of g 2 T are 
important in the integral determining the quantity under consideration. An important 
example for such a quantity is the damping rate of a hard parton, considered in detail 
in the next chapter. So far there are only preliminary or phenomenological attempts 
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to improve this situation. 

In the following chapters the formalism introduced here will be applied to various 
interesting quantities of the QGP. The usefulness and the elegance, but also the 
problems of the effective perturbation theory will be exhibited in these applications. 
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Chapter 4 
Damping rates 



The interaction rate of a particle at finite temperature indicates, how often collisions 
of the particle with the particles of the heat bath take place. It is obtained as a 
convolution of the magnitude of the square of the matrix element responsible for 
the collision with the distribution functions of the external particles. The damping 
rate, defined as the imaginary part of the dispersion relation, follows from the imag- 
inary part of the self energy of the particle, as e.g. the gluon damping rate at zero 
momentum in section 2.5. Through cutting rules (Weldon, 1983) the damping and 
interaction rates are closely related. This connection allows for an easy interpretation 
of the physical process responsible for the damping, which showed us already the 
inconsistency of the gluon damping rate, calculated within naive perturbation theory. 

Owing to long range forces of gauge theories, there are more collision partners 
at large distances. Thus the damping rates in QED and QCD are dominated by 
small momentum transfers, leading in general to infrared divergent results if bare 
propagators are used. (In the case of the gluon damping rate at zero momentum, 
calculated in section 2.5 using bare propagators, a finite, although inconsistent result 
has been obtained, because the damping mechanism relies here on the decay of a 
plasma mode without gluon exchange to lowest order.) These infrared singularities are 
screened by medium effects that can be taken into account using effective propagators 
and vertices in the computation of the matrix elements or the self energies. Hence 
interaction and damping rates are typical examples, requiring the application of the 
Braaten-Pisarski method. 

First we will demonstrate the calculation of the interaction rate of an energetic 
particle (muon, electron, quark, or gluon) in detail. This problem represents the 
simplest application of the Braaten-Pisarski method in gauge theories, since only the 
effective gauge boson propagator is needed. Moreover, these rates are the starting 
point for other interesting quantities, such as the energy loss (chapter 5) or the vis- 
cosity (chapter 6). Furthermore, these rates are already plagued with the problem 
of the missing magnetic screening. Because of these properties damping rates have 
been studied in a large number of publications (Pisarski, 1989b; Lebedev and Smilga, 
1990, 1991, 1992; Thoma and Gyulassy, 1991; Braaten, 1991; Braaten and Thoma, 
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Figure 4.1: Elastic scattering of a muon (bold line) off a thermal electron to lowest 
order. 

1991a; Burgess and Marini, 1992; Rebhan, 1992a; Baier, Kunstatter, and Schiff, 1992; 
Nakkagawa, Niewgawa, and Pire, 1992, 1993; Smilga, 1992; Altherr, Petitgirard, and 
del Rio Gaztelurrutia, 1993; Baier, Nakkagawa, and Niegawa, 1993; Heiselberg and 
Pethick, 1993; Pisarski, 1993; Peigne, Pilon and Schiff, 1993; Thoma, f994a; Baier 
and Kobes, 1994). 

Afterwards we turn to the damping rates of collective plasma modes (Braaten 
and Pisarski, 1990b, 1992b; Pisarski, 1991; Kobes, Kunstatter, and Mak, 1992), such 
as the gluon damping rate at zero momentum, which started the development of the 
Braaten-Pisarski method. The treatment of these rates is much more tedious than 
the one of rates for hard particles, because effective vertices have to be considered 
in addition. An exception is the damping of a soft Yukawa fermion which can be 
calculated easily, as there are no effective vertices in the Yukawa theory (Thoma, 



4.1 Damping rates of hard particles 

Here we will consider the damping rate of a heavy lepton propagating through a rel- 
ativistic electron-positron plasma (Braaten and Thoma, 1991a). This QED problem 
has been investigated in detail, since it exemplifies the application of the Braaten- 
Pisarski method for gauge theories in the simplest possible way. We assume that the 
mass and the momentum of the heavy lepton are much larger than the temperature, 
which in turn is much larger than the electron mass: p, M 3> T 3> m e . This can 
be realized for example for muons at a temperature of about 10 MeV, as it might 
be the case in supernova explosions. We will display the calculation in detail, as it 
reflects the methods used in an ideal manner, and since the rate is closely related to 
observables of the QGP discussed in chapter 5 and 6. 

The interaction rate is caused to lowest order in the coupling constant e by the 
Feynman diagram of fig. 4.1, i.e. by elastic scattering of the muon off the thermal 
electrons and positrons via the exchange of a photon. Here P = (E, p) denotes the 
four momentum of the incoming muon with mass M, energy E = \/p 2 + M 2 , and 
velocity v = p/E. The four momentum K = (fe, k) belongs to the electron in the 
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Figure 4.2: Lowest order contribution to the imaginary part of the self energy of a 
muon (bold line) in naive perturbation theory. 

initial state, where we have neglected the electron mass. Momenta and energies with 
a prime indicate outgoing particles. 

The interaction rate of the process fie" — > fie~ is obtained analogously to the 
zero temperature case, 

x (2tt) 4 6 4 (P + K -P' -K')2(\M\ 2 ) , (4.1) 

where the angle brackets mean averaging over the possible initial and summation over 
the final spin states. The factor of 2 comes from the summation over the spin states of 
the incoming electron. There will be another factor of 2 from the process fie + — > /ie + , 
which has the same matrix element to lowest order. Here we assume an identical 
number of electrons and positrons in the plasma, i.e. a vanishing chemical potential. 
(The generalization of the Braaten-Pisarski method to finite chemical potential has 
been discussed by Vija and Thoma (1994); see also section 7.2.) The only difference 
in the definition of the interaction rate according to (4.1), compared to the zero 
temperature case, is the appearance of the distribution functions. The factor np(k) 
indicates the number of thermal electrons as possible scattering partners, while the 
factor [1 — np{k)] describes Pauli blocking of the scattered electrons. There are no 
distribution functions of the muon which is not a thermal particle. 

Using cutting rules (Weldon, 1983; Keil, 1989) the interaction rate can be obtained 
also from the imaginary part of the muon self energy via 

T(E) = _-L [1 _ n F (E)} tr [(P» lfi + M) Im £(£ + it, p)] . (4.2) 

The self energy corresponding to the matrix element of fig. 4.1 is displayed in fig. 4.2. 
Cutting through the internal fermion lines, shows the equivalence of (4.1) and (4.2). 
(An rigorous generalization of the Cutkosky cutting rules (Itzykson und Zuber, 1980) 
to finite temperature has been given by Kobes and Semenoff (1985, 1986).) 

So far only bare propagators and vertices have been employed (fig. 4.1 and fig. 4. 2). 
The interaction rates, calculated from these diagrams, are proportional to e 4 and 
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Figure 4.3: Lowest order contribution to the imaginary part of the muon (bold line) 
self energy using the effective perturbation theory. 

quadratically infrared divergent. Now we will apply the Braaten-Pisarski method to 
this problem. Intuitively one could replace the bare photon propagator in fig. 4. 1 by 
the effective one. However, the Braaten-Pisarski method has been derived within the 
imaginary time formalism. Therefore we will start from the muon self energy, which 
can be calculated in the imaginary time formalism. 

The lowest order contribution to the imaginary part of the muon self energy, 
using the effective perturbation theory, stems from the diagram in fig. 4. 3. Obviously, 
we do not need an effective muon propagator or vertex in the self energy, since the 
fermion lines are hard due to the muon mass. The effective photon propagator (fig. 2. 6) 
contains an infinite number of electron loops with hard momenta (K ~ T); i.e., the 
damping mechanism is caused by elastic scattering of the muon off the electrons and 
positrons of the plasma via the exchange of a virtual plasma mode. The imaginary 
part of the muon self energy comes from the imaginary part of the effective photon 
propagator, arising for k$ < k 2 . Thus the damping mechanism can be understood as 
virtual Landau damping of the collective photon. 

It is interesting to note that the computation of the muon self energy with an 
effective photon propagator according to fig. 4. 3 is easier than in naive perturbation 
theory, where a two-loop diagram (fig. 4. 2) has to be evaluated. (The one-loop diagram 
with a bare photon propagator contains no imaginary part, since a bare muon cannot 
decay electromagnetically.) 

Here we will follow the derivation of the interaction rate according to Braaten 
and Thoma (1991a). Standard Feynman rules give the following expression for the 
muon self energy of fig. 4. 3: 

v{p) = * J 0j- 4 d iaq) r sin r , (4.3) 

where P' = P — Q . It is convenient to choose the photon propagator in Coulomb 
gauge and to decompose it into longitudinal and transverse contributions as in (2.27). 
Performing the trace over the 7-matrices in (4.2) and employing the imaginary time 
formalism we find 

tr [{P» lfi + M) E*(P)] = 
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go J l Z7r J 

+2 Dj{Q) [p 2 - p q + p • q — (p • q) 2 — M 2 ]} . (4.4) 



The easiest way for evaluating the sum over q is given again by the Saclay method. 
In the case of the muon propagator A(P') we can directly adopt (2.12) and (2.13). 
In the case of the effective photon propagator it is convenient to use the spectral 
representation (Schulz, 1992; Taylor, 1993; Kalashnikov, 1994), 



DIAQ) = f"' ■ (4-5) 

J -co qo — u> + it 

From this the effective photon propagator within the Saclay method can be writ- 
ten as (Pisarski, 1989a): 

pf3 poo 

Di T {q ,q) = - dre qoT dto p L)T {co, q) [1 + n B (to)\ e~" T . (4.6) 

JO J-oo 

(The equivalence of (4.5) and (4.6) follows directly by integrating over r und using 
the expression (2.4) for the Bose distribution.) 

The spectral function contains a contribution from the residue of the pole of the 
effective propagator and a discontinuous contribution from the imaginary part of the 
propagator (Landau damping) if a; 2 < q 2 (Pisarski, 1989a): 

p L , T (u, q) = p r fi(u, q) 6(u - u L , T (q)) + p% T (u, q) % 2 - u 2 ) . (4.7) 

Here LOL^(q) are the poles of the effective propagator, p r £ s T (u,q) the corresponding 
residues, and 

p d ^ = --ImDl T (u } q). (4.8) 

7T 

(The expression for the discontinuous part follows directly by inserting (4.8) into 
(4.5), using the relation Im l/(k — u + ie) = —TrS(k — u), the connection with the 
residue by means of (2.6).) 

In the following we only need the discontinuous part, which is responsible for the 
damping (virtual Landau damping), as discussed above: 

pH'M = ^> L (^)| 2 , 
2q 

\ 3rn 2 u(q 2 -Lo 2 ) 
p T (u,q) = — \D T (u,q)\ . (4.9) 

The summation over q in (4.4) can be done easily using the representation (4.6). 
It gives a factor 8{t — r'), rendering the r' integration trivial. The integration over r 
yields 

d ( P0 ±E>-u)t = ± u .10) 

p ± E> - to y ' 
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Owing to the discrete energies of the muon, p = (2n + l)z7rT, we can set exp(/3p ) = 
— 1 and continue p in the denominators analytically to real values p = E + ie on 
the mass shell afterwards. The imaginary part of the muon self energy is due to the 
imaginary part of the energy denominator 
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So far this expression is exact. Now we will utilize the assumption p } M ^> T . Then 
we can assume the following approximations. First, we may neglect factors such as 
exp(— E/T), exp(— £"/T), np(E) } and riF(E'), which is obvious, since the muon is not 

a thermal particle due to its large mass. Next, we may replace E' = \J(p — q) 2 + M 2 
by E' ~ E — v • q, expanding for small q/E. This expansion is possible because the 
integrand decreases like 1/q 3 ; i.e., only small q contribute to the integral. Then we 
get 6{E — E' — u) ~ 8{u — v • q), which enables us to perform the angle integration: 
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. (4.13) 



Since only small momenta, q <^ T, contribute to the damping rate, also u has to be 
small, to <^ T, allowing for an expansion of the Bose factor: 
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(4.14) 



In the following we only consider the first term of this expansion. The second one does 
not contribute anyway, as the spectral densities are odd functions oiu. The resulting 
integral can be evaluated only numerically. In the the limit M 3> p, however, (4.13) 
can be simplified further, allowing for an analytic evaluation of the rate. From M 3> p 
we get !)<1 and thus u <^ q. Then we may expand the spectral densities for u <^ q 
(Pisarski, 1989b): 
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Substituting these expressions together with (4.14) into (4.13), the interaction rate, 
corresponding to the exchange of a longitudinal photon, is given by 



where we introduced an infrared cut-off, e — > 0, for the ^-integral. 
These results deserve a number of remarks: 

1. The damping rate 7 of a muon, defined as the imaginary part of the dispersion 
relation (Thoma and Gyulassy, 1991), is just one half of the interaction rate, 7 = L/2. 

2. The interaction rate is proportional to e 2 instead of e 4 as expected from naive 
perturbation theory. This anomalous large damping (Lebedev and Smilga, 1990, 
1991, 1992) is caused by the use of the effective photon propagator, similarly to the 
reduction in the coupling constant of the next to leading order term of the effective 
tadpole of the ^-theory, (3.2). (The effective photon masses, m 2 ~ e 2 , appearing 
in the numerator and denominator of the spectral density, cancel each other when 
integrating over q (see also remark 8).) 

3. The logarithmic infrared singularity in the transverse part of the interaction rate 
(4.17) is caused by the absence of a magnetic screening mass; i.e., the rate is sensitive 
to momenta of the order of e 2 T. Since we do not expect a magnetic screening in a 
relativistic electron-positron plasma, this would mean that a muon cannot propagate 
in such a plasma. However, there are doubts (Lebedev and Smilga, 1990, 1991, 1992; 
Smilga, 1992; Thoma, 1994a), whether the damping rate has a physical meaning at 
all. Alternatively the damping could be described by the infrared finite transport 
rate, discussed in chapter 6. On the other hand, a muon in a plasma is never exactly 
on the mass shell due to continuous collisions. If the finite width, i.e. an imaginary 
part of the muon propagator, due to these collisions is taken into account, a finite 
result for the damping rate can be found. Since the width is caused by the damping 
of the muon itself, this method corresponds to a self consistent determination of the 
damping. (Lebedev and Smilga, 1990, 1991, 1992; Altherr, Petitgirard, and del Rio 
Gaztelurrutia, 1993; Pisarski, 1993; Henning, Sollacher, and Weigert, 1994). This 
implies that the infrared cut-off e in (4.17) should be of the order of e. If it is possible 
to derive a consistent result for the damping rate in this way is under debate (Peigne, 
Pilon, and Schiff, 1993; Baier and Kobes, 1994). Another possibility of avoiding the 
infrared singularity of the damping rate by the Bloch-Nordsieck mechanism has been 
proposed by Niegawa (1994). In any case the reduction of the quadratic infrared 
singularity found in naive perturbation theory to a logarithmic one is an important 
progress, allowing for an estimate of the rate up to logarithmic accuracy. 

4. The formula (4.13) also holds for a hard electron with momentum p ~ T, setting 
the velocity equal to v = 1. The calculation, shown above, can be taken over. In 




(4.16) 



The transverse part is logarithmically infrared divergent, 




(4.17) 
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principal there could be a self energy diagram in addition to fig. 4. 3, containing an 
effective electron propagator instead of an effective photon propagator. However, 
owing to momentum conservation at the vertex and the fact that the interaction rate 
is dominated by small momenta of the photon, the electron in the loop of the self 
energy is always hard; i.e., a bare electron propagator suffices. Indeed, it is possible 
to show explicitly that the contribution from the diagram with an effective electron 
and a bare photon propagator is of higher order. The neglect of factors of the form 
exp(—E/T) in (4.13) is approximately possible, as long as the electrons have at least 
thermal energies, i.e. E > 3T. The static approximation, t; <C 1, for the spectral 
densities (4.15), however, has to be revoked. Therefore the integration over u has to 
be performed numerically, yielding Tl = 1.098 e 2 T/47r (Thoma und Gyulassy, 1991; 
Thoma, 1994a); i.e. the deviation from the static approximation is only 10 %. 

5. The interaction rate of a heavy quark, e.g. a charm or bottom quark, in the QGP 
simply follows by replacing the square of the electric charge e 2 of the muon in (4.16) 
and (4.17) by the chromoelectric charge Cf§ 2 of the quark, where Cf = 4/3 denotes 
the Casimir invariant of the fundamental representation of SU(3). In the case of a 
heavy quark at rest, v = 0, this result has been derived by Pisarski (1989b). 

6. The interaction rate of a light quark with hard momentum has been considered by 
means of the Braaten-Pisarski method first by Thoma und Gyulassy (1991). Again 
the electron charge in the rate of a hard electron has to be replaced by the color 
charge of the quark resulting in T ~ (Cpg^T /2tt) ln(l/g) to logarithmic accuracy, 
where an infrared cut-off ^ < 1 has been chosen. In a similar way this result (up to 
a factor of 3) was found already earlier by Lebedev and Smilga (1990, 1991), which 
was verified later on also by Burgess and Marini (1992) and Rebhan (1992a). 

7. The interaction rate of a hard gluon follows from the one of a light quark by 
replacing Cf by the Casimir invariant of the adjoint representation, Ca = 3. This can 
be seen either by explicit evaluation of the gluon self energy using an effective gluon 
propagator (Braaten, 1991; Burgess und Marini, 1992) or immediately by comparing 
the scattering diagrams of elastic quark scattering (fig.4.1) with the ones of elastic 
gluon scattering, caused by the scattering of the gluon off the thermal quarks and 
gluons via the exchange of a gluon. 

8. The interaction rates of light partons are independent of the parton momentum p 
in the limit p ^> lo } q. Moreover, they do not depend on the number of thermalized 
flavors in the QGP, since the effective gluon mass m g in the spectral density drops 
out after integrating over q. Though the number of possible scattering partners is 
increased if additional quark flavors are present, at the same time the screening is 
enlarged. Both effects cancel each other exactly (Thoma, 1994a). 

9. The quark damping rate has been calculated also in covariant gauges (Baier, Kun- 
statter, and Schiff, 1992). The gauge dependent part proportional to the gauge fixing 
parameter £ is infrared divergent and has to be removed by an infrared regulator. 
The gauge dependent part vanishes only if one first puts the self energy on mass 
shell, before the regulator is set equal to zero. Interchanging both the limits, a gauge 
dependent contribution survives. This observation caused doubts about the gauge 
independence of the Braaten-Pisarski method, which, however, could be disproved 
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meanwhile by arguing for the correct order of the limits (Rebhan, 1992b; Braaten 
and Pisarski, 1992b; Nakkagawa, Niegawa, and Pire, 1992, 1993; Altherr, Petitgirard, 
and del Rio Gaztelurrutia, 1993; Baier, Nakkagawa, and Niegawa, 1993; Kobes and 
Mak, 1993). 

10. In a QGP there might be the possibility of a magnetic screening mass due to 
non-perturbative field configurations, e.g. chromomagnetic monopoles. This mass 
should be proportional to g 2 T (Linde, 1980; Gross, Pisarski, and Yaffe, 1981) and 
could serve as an infrared regulator for the interaction rate (Pisarski, 1989b). Indeed, 
a magnetic gluon mass has been found in lattice calculations (Biloire, Lazarides, and 
Shafi, 1981; deGrand and Toussaint, 1982) as well as in classical considerations (Biro 
and Miiller, 1993). Despite of the completely different nature of both the approaches, 
identical results m mag = 0.3lg 2 T have been obtained. Pisarski (1993) combined the 
magnetic mass with the finite width of the quark phenomenologically. Extrapolating 
these results to realistic values of the strong coupling constant, e.g. a s = 0.3, gives 
T q ~ 0.3T for quarks and T g ~ 0.7T for gluons (Thoma, 1994a). 

11. Assuming that the interaction rates of light and hard partons describes typical 
relaxation and thermalization times of the QGP, r gj3 = l/r ?)S , surprisingly short 
thermalization times of about T q ~ 1 - 2 fm/c und T g ~ 0.5 - 1 fm/c are obtained 
choosing a temperature of 300 MeV (Thoma, 1994a). These results agree with the 
ones found in numerical simulations of ultrarelativistic heavy ion collisions (Wang and 
Gyulassy, 1991; Geiger and Miiller, 1992; Geiger, 1993). Furthermore one observes 
that the gluon component of the parton gas thermalizes faster than the quark compo- 
nent (Shuryak, 1992; Geiger, 1993; Biro et al., 1993) due to the stronger interactions 
of gluons compared to quarks. Whether the interaction rates, however, determine the 
thermalization times directly is questionable, as discussed in chapter 6. 

12. The interaction rate determines the color conductivity and diffusion of the QGP, 
which describes the relaxation from a deviation of the color equilibrium distribu- 
tion (Selikhov and Gyulassy, 1993; Heiselberg, 1994b). The color conductivity, for 
instance, is given by o c = 2m 2 , /T g (Selikhov and Gyulassy, 1993). 

13. The rates of inelastic processes such as gg — > qq or gg — > ggg have not been 
treated by means of the Braaten-Pisarski method up to now. Because they are of 
interest for the chemical equilibration of the QGP (Biro et al., 1993) and the radiative 
energy loss of energetic partons (Gyulassy and Wang, 1994), they have been estimated 
using bare Green's functions and the effective gluon mass m g and quark mass m q as 
infrared cut-offs (Xiong and Shuryak, 1992; Biro et al., 1993; Altherr and Seibert, 
1994). In this way the electric and magnetic infrared singularities are removed at the 
same time, of course without any justification. 

14. The mechanism of photon damping to lowest order is caused by Compton scat- 
tering as opposed to elastic scattering, which is responsible for parton and lepton 
damping. The photon damping rate will be discussed in detail in chapter 7. 
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Figure 4.4: Self energy contributions to the damping of a soft gluon. 

4.2 Damping rates of soft particles 

Now we will turn to the damping rates of particles with soft momenta, i.e. of the 
order </ /'. Such particles can be regarded as collective modes, e.g. plasmons, as 
discussed in section 2.4. A consistent treatment of those rates again requires the use 
of the Braaten-Pisarski method. Again the rates follow from one-loop self energies 
containing effective Green's functions. Since the external momenta are soft now, all 
lines of the diagram determining the self energy can be soft. Thus effective vertices 
as well as propagators have to be used. 

4.2.1 Gluon damping rate 

The leading order contribution to the gluon damping rate at zero momentum 
stems from the diagrams of the gluon self energy in fig. 4. 4. Owing to the complicated 
expressions for the effective vertices and propagators, the evaluation of these diagrams 
is very tedious. Since we only need the imaginary part of these diagrams at zero 
momentum (see (2.31)), the calculation can be simplified considerably. After all, the 
calculation is far from being trivial. The gauge independent result of this lengthy 
calculation (Braaten und Pisarski, 1990b), 



represents the solution of the plasmon puzzle discussed in section 2.5. The damping 
rate is significantly larger than the one found by using naive perturbation theory in 
temporal or Coulomb gauge, (2.34) with a = 1, and positive in contrast to the result, 
obtained using bare Green's functions in covariant gauge. Hence there is no indication 
of an instability of the QGP. 

The physical damping mechanism can be recognized again by cutting the diagrams 
in fig. 4. 4. There are the following possibilities: We can cut through two lines ("pole- 
pole term"), through one line and one blob ("pole-cut term") or through two blobs 
("cut-cut term"). Cutting through a blob means cutting a HTL self energy, leading 
to virtual Landau damping. The pole-pole term corresponds to the decay of the 
collective mode (at zero momentum) in two collective modes, which is forbidden 
by kinematics. The pole-cut term describes gg —* gg scattering and the cut-cut 
term inelastic gg — > ggg processes (bremsstrahlung). So the damping mechanism of 
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Figure 4.5: Self energy contributions to the damping of a soft quark. 

collective gluons at zero momentum relies on elastic and inelastic scattering processes 
between plasma waves and thermal gluons, as opposed to the decay of the collective 
mode in bare gluons as expected in naive perturbation theory (see section 2.5). This 
observation demonstrates once more that the Braaten-Pisarski method leads not only 
to quantitatively improved results, compared to naive perturbation theory, but also 
contributes to a deeper insight into the physical processes in the QGP. 

It is interesting to note that, in contrast to the hard rates, there are no infrared 
singularities due to the missing static magnetic screening. Although the rates at zero 
momentum are caused by scattering, the exchange of a transverse gluon does not 
lead to a divergence, since magnetic processes are suppressed for vanishing external 
momentum. 

Finally, we would like to note that investigations of the chaoticity of a classical 
Yang-Mills system on the lattice, using the Hamilton formalism in temporal gauge, 
yielded a Lyapunov exponent which surprisingly agrees with the gluon damping rate 
at zero momentum, (4.18), (Miiller and Trayanov, 1992; Biro et al., 1994). A possible 
connection between the damping rate and the Lyapunov exponent might rely on the 
fact that both quantities measure the entropy increase of the system (Gong, 1993). 

4-2.2 Quark damping rate 

The quark damping rate at zero momentum follows from (Braaten and Pisarski, 
1992b) 

7±(0) = -~ tr[ la /mS>, + it. 0)] , (4.19) 

o 

where j± denotes the rates for positive and negative helicity eigenstates, respectively, 
which are identical at zero momentum analogously to the transverse and longitudinal 
gluon damping rates, and m q = gT / 1/6 is the thermal quark mass. For the quark self 
energy E* we have to adopt the diagrams in fig. 4. 5. Here an effective vertex between 
two quarks and two gluons which has no bare counterpart is needed. Independently 
Kobes, Kunstatter, and Mak (1992) and Braaten and Pisarski (1992b) found 

7±(0) = &f^, (4.20) 

1Z7T 

where the constants b = 5.63 and 5.71 apply to two and three thermal quark flavors, 
respectively. The damping mechanism is caused by elastic (2 — > 2) and inelastic 
(2 —> 3) scattering processes, as in the case of the gluon damping. 
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Chapter 5 
Energy loss 



The transition and energy loss of a charged particle through matter has been investi- 
gated intensively since decades (Jackson, 1975). For example the famous Bethe-Bloch 
formula describes the energy loss of a charged particle with a mass, which is large 
compared with the electron mass, in non-relativistic matter, such as a plasma, by 
collisions of the particle with the electrons of the matter. In the case of relativistic 
particles radiative corrections, i.e. energy loss by bremsstrahlung, become important. 
This process dominates especially for light particles such as electrons. 

5.1 Motivations and estimates 

The energy loss of particles in a relativistic medium, e.g. a QGP or a relativis- 
tic electron-positron plasma, requires a quantum field theoretic treatment. For e.g. 
cross sections of quark scattering, following from QCD, enter the energy loss. In this 
way Bjorken (1982) estimated the energy loss of energetic partons in the QGP. This 
quantity is of interest, as it is closely related to the appearance of jets in ultrarela- 
tivistic heavy ion collisions, which are supposed to be observable at RHIC and LHC 
for the first time in heavy ion collisions (Eskola, Kajantie, and Lindfors, 1989). 

Jets can be considered as an external probe for the fireball. For in the case of 
a relativistic heavy ion collision the energetic parton, produced in primary parton 
collisions has to traverse the fireball (fig. 5.1), being decelerated and deflected. This 
modification by the medium leads to an asymmetric jet distribution, while in electron- 
positron or proton-proton collisions jets originating from a quark-antiquark pair are 
emitted in opposite directions. Indeed, the deflection of the jet axis (acomplanarity) 
has been discussed as a possible signature of the QGP formation (Apple, 1986; Blaizot 
and McLerran, 1986; Rammerstorfer and Heinz, 1990). Unfortunately, this deflection 
may also arise in a hadronic fireball, rendering the distinction of the two phases by 
means of the jet-acomplanarity impossible (Rammerstorfer and Heinz, 1990). 

Rather the deceleration of jets might serve as a signal for the QGP. According 
to preliminary estimates the energy loss of energetic partons in the two phases might 
differ by an order of magnitude; i.e., jets in hadronic matter might be suppressed as 
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Figure 5.1: Jets in ultrarelativistic heavy ion collisions. 

opposed to jets in a QGP ("jet unquenching" ) (Gyulassy and Pliimer, 1990). Whereas 
the energy loss of quarks in hadronic matter is given by the string tension of about 
1 GeV/fm, Bjorken's estimate (Bjorken, 1982) for a quark with an energy of about 
E = 10 - 20 GeV in the QGP yielded dE/dx = 0.1 - 0.2 GeV/fm. 

Bjorken's estimate is based on the use of elastic qq- and gg-cross sections to lowest 
order perturbation theory, i.e. on diagrams of the type shown in hg.4.1. The infrared 
singularities in these cross sections have been removed similarly as in the case of 
the inelastic rates, discussed in the last chapter, by introducing the Debye mass, 
m 2 D = 3m^, as infrared regulator. 

However, according to Fermi (Jackson, 1975; Ichimaru, 1973) the infrared contri- 
bution to the energy loss due to large distance collisions can be described as follows: 
The incoming charged particle induces an electric held in the medium, which reacts 
via the Lorentz force upon the incoming particle causing its deceleration. By means 
of this classical consideration the energy loss can be related to the dielectric constant 
of the medium. In the case of quarks one finds (Thoma and Gyulassy, 1991) 

§=-£&Jt S-l ^ u i /m EZ ' + ( " 2 e ~ " 2) Im ^ tT ~ ,2r '] • (5 - 1) 

where Cp = 4/3 and v is the velocity of the quark. Adopting the expressions (2.18) for 
the longitudinal and transverse dielectric constant together with the HTL gluon self 
energy (2.16), a finite result has been found from (5.1), rendering the introduction of 
an infrared cut-off obsolete. Mrowczyhski (1991) combined the infrared contribution 
found in this way with the one caused by short distance collisions, estimated by 
Bjorken, introducing the Debye mass as a separation scale. A similar proceeding has 
been proposed by Koike and Matsui (1992). In both cases the final result depends on 
the somewhat arbitrary separation scale. 

The energy loss of a quark in the QGP, following from these considerations, can 
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transfer in the parton collision. Assuming a s = 0.2, Nj = 2, T = 250 MeV, and 
E = 20 GeV, as chosen by Gyulassy and Pliimer (1990), an energy loss of about 0.2 
GeV/fm results. The small value a s = 0.2, motivated by lattice calculations (Gao, 
1990) and renormalization scheme arguments (Huang and Lissia, 1994), however, is 
affected by a large uncertainty. Already for a s = 0.4 the energy loss according to (5.2) 
is of the order of the string tension. Also the choice of q ma x and q m in is ambiguous. 
Assuming e.g. q ma x = E results in a doubling of the energy loss. Finally, another 
important uncertainty is due to the neglect of higher order effects, such as radiation. 

The uncertainty in q m in can be removed employing the Braaten-Pisarski method, 
the one in q ma x taking into account the kinematics rigorously at least in the case of 
heavy quarks. In the forthcoming sections the energy loss of heavy quarks due to 
elastic scattering in the QGP will be derived consistently, leading to a result which is 
infrared finite, gauge independent, complete to leading order in the coupling constant, 
and independent of arbitrary parameters like cut-offs and separation scales. 

For estimating the energy loss we also will assume a s = 0.2, as there is no reliable 
method for determining the temperature dependent effective coupling constant up 
to now. (Attempts for deriving a temperature dependent running coupling constant 
starting from the renormalization group equation at finite temperature (Baier, Pire 
and Schiff, 1991; van Eijck, Stephens and van Weert, 1994) did not lead to unam- 
biguous results so far.) The validity and limitation of the extrapolation of the results 
for dE/dx } obtained by means of the Braaten-Pisarski method for g <C 1, to realistic 
values of g = 1.5 - 2.5 will be discussed in section 5.5. 

The calculation of the energy loss by bremsstrahlung is much more involved, since 
higher order diagrams as well as new medium effects (Landau-Pomeranshuk effect) 
arise. Consequently this contribution to the total energy loss is not yet treated within 
the Braaten-Pisarski method. In section 5.6 rough estimates will be given. 

Besides for jet quenching as a possible signature for the QGP, the energy loss 
has been used as an input for the numerical simulation of nucleus-nucleus collisions 
(Wang and Gyulassy, 1992), for estimating thermalization times (Shuryak, 1992) in 
ultrarelativistic heavy ion collisions, and for the mean free path of heavy quarks in a 
primordial QGP (Dine et al., 1992). Another quantity, closely related, is the diffusion 
rate of a charm quark, which plays a role for the J/^-suppression. It has has been 
calculated by Svetitsky (1988) analogously to the energy loss according to Bjorken. 

5.2 Quantum field theoretic definition 

For simplicity sake we will first consider the energy loss of a muon in a relativistic QED 
plasma. Furthermore we will restrict ourselves to the collisional energy loss, which 
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is given to lowest order by elastic scattering of the muon off the thermal electrons 
and photons of the plasma via the exchange of a photon (Fig.4.1). (Although elastic 
scattering does not cause an energy transfer in the center of mass system, it does in 
the system of the heat bath.) 

A quantum field theoretic definition of the energy loss can be obtained in the 
following way (Braaten and Thoma, 1991a): The mean distance between two collisions 
of a muon with velocity v in the plasma is given by Ax = v/T } where T is the 
interaction rate of section 4.1. (As usual, we assume that the dimension of the plasma 
is large compared to the mean free path of the muon, corresponding to a large number 
of collisions, and that the energy transfer per collision is small compared to the energy 
of the muon.) The energy loss of a muon per collision is given by AE = E — E' = lo } 
where u denotes the energy transfer by the exchanged photon with four momentum 
Q = (cu,q). The mean energy loss is now obtained by averaging over the interaction 
rate Y times the energy transfer u and dividing by the velocity v of the muon. This 
is expressed in a symbolic formula as 

^ = - [dTuj, (5.3) 
ax v J 

meaning that the photon energy u has to be inserted under the integrals that define 
the interaction rate, i.e. (4.1) or (4.2). The energy loss can be calculated analogously 
to the interaction rate either from the scattering matrix element or from the imaginary 
part of the muon self energy. 

However, there is an essential difference compared to the interaction rate: The ad- 
ditional factor to in the energy loss changes the infrared behavior completely. Whereas 
the interaction rate is quadratically infrared divergent using a bare photon propaga- 
tor and logarithmically using a resummed one, the energy loss is logarithmically in- 
frared divergent in naive perturbation theory but finite applying the Braaten-Pisarski 
method. This can be seen by inspecting (4.13). If u is inserted under the integral 
there, the infrared divergence is reduced directly by the factor u and in addition by 
the fact that now the second term of the expansion of the distribution function in 
(4.14) contributes, because the spectral functions pf; s T (ijj,q) are odd functions of to. 
Indeed, the use of the effective photon propagator suffices for producing a finite result, 
although the effective photon propagator contains no static magnetic screening. This 
dynamical screening of infrared singularities for quantities that are logarithmically 
divergent in naive perturbation theory has been observed first in the case of the vis- 
cosity (Baym, Monien, and Pethick, 1988; Pethick, Baym, and Monien, 1989; Baym 
et al., 1990). 

Owing to the additional factor u, a further modification arises compared with 
the interaction rate besides the different infrared behavior. While in the latter case 
only soft momentum transfers, g~eT, need to be considered because of the strong 
decrease of the integrand with momentum, T ~ / dq/q 3 } all kinematically possible 
momentum transfers have to be taken into account in the case of the energy loss, 
dE I dx ~ / dq/q. Since the effective photon propagator in the muon self energy, (4.3) 
and fig. 4. 3, contributes to damping only for soft momenta - for hard momenta it 
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reduces to the bare propagator, which contains no imaginary part - it is not possible 
to compute the energy loss by simply using the effective propagator in the muon self 
energy for the entire momentum range. 

An alternative would be the use of the effective photon propagator in the scatter- 
ing matrix element according to fig. 5. 2. This diagram reduces to fig.4.1, as desired, for 
large momentum transfers Q and considers screening effects for soft momenta at the 
same time. This approach, however, suffers from two problems: First, it cannot be 
derived rigorously from the Braaten-Pisarski method, which relies on the imaginary 
time formalism; i.e., strictly speaking, the use of effective propagators is only possible 
in loop diagrams. (However, both approaches give the same result in the case of the 
energy loss for soft momentum transfers (Braaten and Thoma, 1991a).) Secondly, 
the diagram in fig. 5. 2 integrated over the external momenta contains higher order 
contributions in the coupling constant, as can be seen, studying a simple analytic 
example (Thoma, 1994d): 



Here t = Q 2 stands for the momentum of the photon with maximum momentum 
transfer s = (P + K) 2 and A* = l/(t + Ji 2 ) for the effective propagator with thermal 
mass Jj ~ eT. The factor t under the integral softens the infrared divergence, appear- 
ing for Ji, = 0, as in the case of the energy loss. The second term in the parentheses of 
(5.4) yields contributions to higher orders, 0(e 6 ). In the simple example above these 
contributions can be removed by neglecting these terms making use of s 3> Ji 2 . How- 
ever, if an integration over distribution functions (see e.g. (4.1)) has to be performed, 
an analytic result as in (5.4) cannot be given. Hence the higher order contributions 
cannot be separated off easily. Consequently, the final result is not consistent in the 
order of the coupling constant. If we want to obtain a complete expansion up to order 
e 6 , we would have to consider additional complicated diagrams, e.g. with effective 
vertices. 

5.3 Braaten-Yuan prescription 

A method for the consistent treatment to leading order of quantities which are loga- 
rithmically infrared divergent, using bare propagators, has been proposed by Braaten 
and Yuan (1991). Introducing a separation scale q*, similarly as in the derivation 
of the photon self energy in the HTL approximation in appendix B, restricted by 
eT <C q* <C T but otherwise arbitrary, soft and hard momentum transfers can be 
treated separately. 

For soft momentum transfers, q < q* } we adopt the Braaten-Pisarski method for 
example in the case of the energy loss by starting from the muon self energy with an 
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Figure 5.2: Elastic scattering of a muon (bold line) off a thermal electron via the 
exchange of a collective photon. 

effective photon propagator (fig. 4. 3). Neglecting terms of the order of eT compared 
to q*, wherever it is possible, i.e. causing no divergence, we obtain an expression of 
the form B lxi(q*/eT) + A so j t with constants A so j t and B that are independent of q*. 
The factor of eT under the logarithm is due to the use of the effective propagator. 

The hard contribution (q > q*) follows from using only bare Green's functions for 
example in the matrix element and from taking account of q* «C T. In this way an 
expression of the form B ln(T/g*) + Ah ar( i is achieved, where B is the same constant 
as in the soft part. Adding up both the contributions, B ln(I/e) + A so j t + Ah ar( i is 
found; i.e., the arbitrary separation scale drops out, reflecting the consistency of the 
final result to leading order. 

Applying this prescription to our simple example (5.4) we get I so ft = e 4 [ln(t* / p 2 ) — 
1] if t < t* and Ihard = e 4 ln(s/t*) if t > t* using the bare propagator A(t) = 1/t. 
In total / = e 4 [ln(s/p 2 ) — I] follows, which agrees with (5.4) in the limit s 3> Jj 2 
and is not contaminated by higher orders (Thoma, I994d). The advantage of the 
Braaten-Yuan prescription is the possibility to separate off the higher order terms 
analytically, even in cases where numerical integrations over distribution functions 
are required, as opposed to using the effective propagator in the matrix element for 
all momentum transfers. In conclusion, the Braaten-Yuan prescription provides an 
analytic way for obtaining results that are complete to leading order in the coupling 
constant. 



5.4 Energy loss of muons 

First, we will consider the soft part, q < q* } combining (5.3) with (4.2) and (4.3). 
Proceeding analogously to the interaction rate but using the second term in the 
expansion of the Bose enhancement factor in (4.14), since the spectral functions are 
odd functions of tu, we find 



dx 
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dq q I did to 
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dis 



. (5.5) 
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Figure 5.3: Compton scattering from a muon (bold line). 

The expression (5.5) coincides with the one, (5.1), obtained by Thoma and Gyu- 
lassy (1991), if (5.1) is carried over to the QED case by substituting CfQ 2 by e 2 , and 
if (2.18), (2.26), and (4.8) are combined, yielding (Thoma and Gyulassy, 1991) 

j — 1 2 dis 

Ime L = -irq p L , 
Imi^tT-q 2 )- 1 = -7T pf . (5.6) 

This agreement demonstrates the equivalence of the HTL approximation, which is 
identical to the high temperature limit, with the classical approximation, used in 
deriving (5.1). 

Changing variables from q and u to q and x = co/q } the integration over q can be 
done exactly. The ^-integral, however, has to be evaluated numerically. Considering 
furthermore q* ^> eT, 

e4 T 2 f 1 l~v\ l + v\ fa* . \ 

In- h n JL. + A S0/i (t)) (5.7) 




24tt \v 2v 2 1 - v / V eT 

sojt \ / 

is obtained, where A so j t (v) begins at A so j t (0) = 0.049, increases to a maximum of 
^4 so /t(0.96) = 0.292 monotonically and decreases then monotonically to A so j t (l) = 
0.256. 

As already mentioned above, the soft part of the energy loss can also be derived 
from the matrix element with an effective photon propagator according to fig. 5. 2 
(Braaten and Thoma, 1991a). This observation confirms the conjecture that the 
imaginary part of the diagram of fig. 4. 3 corresponds to elastic scattering via the 
exchange of a collective plasma mode, meaning the extension of the cutting rules 
(Weldon, 1983) to collective modes, not proven in general so far. 

The hard part, q > q* , can be achieved starting from the matrix element combin- 
ing (5.3) and (4.1) and using bare Green's functions. Besides electron-muon scatter- 
ing, the Compton processes shown in fig. 5. 3 contribute to the same order in e. These 
contributions, however, are suppressed by a factor (T/M) 2 (Braaten and Thoma, 
1991a), where M is the mass of the muon, assumed to be large compared to the 
temperature. 

The matrix element for elastic scattering of muons off electrons and positrons, 
averaged over the initial states and summed over the final ones, reads (Braaten and 
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Thoma, 1991a) 



e 4 

(\M | 2 ) = 8 — (P ■ K P' ■ K' + P ■ K' P' ■ K - M 2 K ■ K') . (5.8) 

This result has been obtained in Feynman gauge, whereas the soft part was calculated 
using Coulomb gauge, which causes no problem because of the gauge independence 
of the hard as well as the soft part in our treatment. 

Inserting (5.8) in (5.3) together with (4.1), a somewhat lengthy calculation gives 
(Braaten and Thoma, 1991a) 

e 4 T 2 (I 1-v 2 , l + v\ / T E A . \ , , 

ln- + ln— + A hard {v)) , (5.9) 

hard 




24tt \v 2v 2 1-vl \ q* M 



where we assumed p } M ^> T again. The maximum momentum transfer, providing 
the upper limit for the integration in deriving (5.9), has been determined from the 
kinematics of the scattering, containing no arbitrary assumption, in contrast to the 
one used by Bjorken in the case of light partons. The function Ahardiy) decreases 
monotonically from Ah ar( i(0) = 1.239 to Ah ar( i(l) = 1.072. 
Combining (5.7) and (5.9), we end up with 

dE e 4 T 2 (I 1-v 2 , l + v\ / E 1 \ 

In— ^ In— + ln- + A(v)) . (5.10) 



dx 247r \v 2v 2 1 - v J V M e 

As expected according to the Braaten-Pisarski method, the final result indeed is 
independent of the separation scale q*. The function A(y) = A so j t (v) + Ahardiy) 
increases monotonically from A(0) = 1.288 to A(0.88) = 1.478 and decreases then 
monotonically to A(l) = 1.328. In fig. 5. 4 the result (5.10) is compared with the 
ones by Bjorken (1982) (dotted line) and Thoma and Gyulassy (1991) (dashed line), 
carried over to the energy loss of a muon. In contrast to the complete result (5.10) the 
estimates of Bjorken and Thoma and Gyulassy, respectively, depend logarithmically 
on M/T, for which M/T = 10 has been chosen in fig. 5. 4. 

The energy loss (5.10) is proportional to e 4 instead of e 2 , as it is the case for the 
damping rate due to the different infrared behavior. The factor of e 2 in the effective 
photon mass entering the denominator of the effective propagator shows up under the 
logarithm in the case of the energy loss, whereas it reduces the order in the coupling 
constant from e 4 to e 2 in the case of the damping rate. 

In the limit of small and large velocities (u — > and v — > 1, respectively) the 
formula (5.10) is not valid anymore, since the approximations v 3> T/E and E <^ 
M 2 /T } respectively, assumed in deriving (5.10) do not hold anymore. After all it is 
possible to derive analytic expressions analogously to (5.10) in both these cases. For 

' dE ^ ^ 'loI + J.288), (5.11) 



dx „ 12irMv \ e 
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Figure 5.4: Energy loss of a muon in dimensionless units as function of the muon 
velocity according to (5.10) (solid line), to Thoma and Gyulassy (1991) (dashed line), 
and to Bjorken (1982) (dotted line). 

has been found (Braaten and Thoma, 1991a), where the negative sign means that 
muons with subthermal energies gain energy from the plasma. 
For v — > 1 Braaten and Thoma (1991a) obtained 



Finally the expression (5.10) describing the energy loss of a massive charged 
particle in a relativistic plasma shall be compared with the one in a non-relativistic 
plasma (Bethe-Bloch-formula), which reads (Jackson, 1975) 



Here u v i denotes the plasma frequency (thermal photon mass) and m e the electron 
mass, and A a constant of the order of one. The non-relativistic formula (5.13) shows 
a minimum as a function of the velocity, whereas the relativistic one, (5.10), increases 
monotonically with v. 

5.5 Energy loss of partons 

The energy loss of a heavy quark Q with mass M 3> T and momentum p ^> T, e.g. 
of a charm or bottom quark, in a QGP can be treated analogously to the energy loss 






(5.13) 
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Figure 5.5: Scattering of a heavy quark (bold line) off thermal quarks and gluons to 
lowest order. 



of a muon (Braaten and Thoma, 1991b) . The processes responsible for the collisional 
energy loss are scattering processes of the heavy quark off the thermal quarks q and 
gluons g, displayed in fig. 5. 5. In contrast to the QED case, Compton scattering cannot 
be neglected, since these contributions (s- and u-channel) do not cancel anymore in 
the limit T/M — > due to color factors. In addition, there is a new diagram describing 
the elastic scattering of the quarks off a gluon (t-channel). 

Analytic results can be found, as in the case of the muon, for the energy regimes 
<C M 2 /T and E 3> Af 2 /T, for which the maximum momentum transfer has a 
simple form, allowing for an analytic integration in the hard part. 



First we will consider the limit E «C M /T. The soft contribution (q < q 



is 



obtained from the one of the muon (5.7) by replacing e by g } the thermal photon mass 
m 7 = eT/3 by the thermal gluon mass m g = (^Jl + Nf/6/\/3) //'/'. and multiplying 
with a color factor Cp = 4/3. These modifications follow immediately comparing the 
quark self energy diagram, responsible for the soft part (fig. 5. 6), with the one of the 
muon self energy (fig. 4. 3) using the standard Feynman rules for QED and QCD. 

The hard part (q > </ ' ) due to Qg-scattering in fig. 5. 5 follows from the hard part 
of the muon energy loss by replacing e by </. multiplying with a color factor of 2/3, 
and summing over the possible Nj flavor states of the incoming thermal quark. The 
hard contribution coming from Qg-scattering has to be calculated from the corre- 
sponding diagrams in fig. 5. 5. Thereby we may set q* equal to zero in the Compton 
contributions, since the s- and u-channels do not lead to infrared singularities. The 
explicit expressions for the hard part, caused by the diagrams of fig. 5. 5, are presented 
in the literature (Braaten and Thoma, 1991b). 

Combining the soft contribution with the hard one, q* cancels again, as expected, 
and the final result for the energy loss of a heavy quark for E <C M 2 /T is given by 

dE 8ira 2 T 2 ( N f \ (l 1 - v 2 , 1 + v\ f A , s ET \ 

'' 1,1 ln^^ In [2 N f +6 B{v) — , (5.14) 



dx 3 V 6 / \v 2v 2 1 — v J \ m g M 

where B(v) is a smooth function of the quark velocity v that increases monotonically 
from B(0) = 0.604 to i?(0.88) = 0.731 and decreases then monotonically to B(l) = 
0.629. 

The energy loss in the limit E ^> M 2 /T can also be carried over partly from the 
one of a muon for v 1, (5.12). Solely the Qg-contribution has to be calculated 
explicitly, where now Compton scattering can be neglected. In this way the energy 
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Figure 5.6: Self energy contribution to the soft part of the energy loss of a heavy 



loss of a heavy quark with E 3>> M 2 /T has been derived (Braaten and Thoma, 1991b), 



Demanding a continuous transition from (5.14) to (5.15), a crossing energy E cross = 
1.80Af 2 /T, if Nj = 2, follows. Applying these results to realistic situations, we choose 
a typical temperature of T = 250 MeV and extrapolate the results (5.14) and (5.15), 
derived in the weak coupling limit, g <C 1, to a s = 0.2. The energy loss of a charm 
quark ( .\ I = 1.5 GeV) and a bottom quark (M& = 5 GeV) are plotted in hg.5.7 and 
fig. 5.8, respectively, as a function of the quark momentum. The crossing energy in 
the case of a charm quark is E cross = 16 GeV leading to the kink in hg.5.7. Since 
this kink is not very pronounced, the use of (5.14) below E cross and of (5.15) above 
it appears to be a good approximation. In the case of a bottom quark, for which 
Ecross = 180 GeV, only the expression for E <C M 2 /T is shown. Comparing the 
complete results (5.14) and (5.15) with the estimates of Bjorken (1982) (dotted line) 
and Thoma and Gyulassy (1991) (dashed line), one observes that the hrst one is a 
better approximation for large momenta, while the latter one for small momenta. In 
the case of a charm quark with a typical energy of E = 20 GeV the energy loss is 
given by about dE / dx ~ 0.3 GeV/fm, whereas the energy loss of a bottom quark is 
just about half of it. This suggests that in particular jets with a bottom quark as 
leading particle are almost unquenched in a QGP. 

Finally, we will point out a peculiarity of the curves in hg.5.7 and fig. 5. 8. For 
small momenta, p < 1.2 GeV for the charm quark and p < 4.1 GeV for the bottom 
quark, the energy loss becomes negative. This cannot be explained as an energy 
gain at subthermal energies, which are considerably smaller. The reason for this 
unphysical result relies on the daring extrapolation from g <C 1 to g = 1.6. As a 
matter of fact, the formula (5.14) fails if // > 1.08 and E < 0.93gM hold at the 
same time, since the logarithm in (5.14) then becomes negative. This behavior is not 
immediately obvious, as the energy loss is derived from the magnitude of the square of 
a matrix element. (This statement also applies to the soft part, since it can be derived 
from the matrix element with an effective propagator as well.) However, using the 
Braaten- Yuan prescription, positive terms have been omitted, namely terms which 
do not diverge for q* — > oo in the soft part and for q* — > in the hard part. These 
terms can be neglected if g <C 1 and, actually, have to be neglected in order to derive 
a consistent final result, as discussed in sections 5.2 and 5.3. Here a limitation of the 



quark (bold line). 




(5.15) 



48 



> 

CD 

CD 

x 

T5 

LU 

T3 




p [GeV] 

Figure 5.7: Energy loss of a charm quark as function of the quark momentum accord- 
ing to (5.14) and (5.15) (solid line), Thoma and Gyulassy (1991) (dashed line), and 
Bjorken (1982) (dotted line). 




p [GeV] 

Figure 5.8: Energy loss of a bottom quark as function of the quark momentum ac- 
cording to (5.14) (solid line), Thoma and Gyulassy (1991) (dashed line), and Bjorken 
(1982) (dotted line). 
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application of the Braaten-Pisarski method to leading order in g becomes apparent. 
For large momenta, however, p ^> T a reasonable behavior is found, which indicates a 
justification for the application of the Braaten-Pisarski method to realistic situations. 

Next we consider the energy loss of a light quark or gluon (Thoma, 1991b). Simi- 
larly to the energy loss of an electron in non-relativistic matter it is more complicated, 
because light partons may lose their entire energy in a single collision, causing large 
fluctuations in the energy loss. Moreover, new gg-scattering diagrams (u-channel), 
leading to infrared singularities for large momentum transfers, arise. Furthermore, the 
integrations over the matrix elements in the hard part cannot be performed analyti- 
cally, since the assumption M 3> T, used for deriving (5.9), does not hold anymore. 
On the other hand, momentum transfers of the order of the momentum of the in- 
coming particle do not contribute to jet-quenching, since then the energy and the 
momentum of the incoming parton will be transferred to a thermal parton, taking 
over the role of the leading particle in the jet. Consequently, the energy stays within 
the jet and will not be dissipated into the QGP. In order to take this behavior into 
account, we introduce a maximum momentum transfer q, which is given by about 
q = p/2; i.e., we neglect ultrahard momentum transfers for describing the physically 
interesting energy loss of a jet in the QGP. Then it can be shown (Thoma, 1991b) 
that this energy loss follows from the one of a heavy quark in the E 3> M 2 /T limit 
simply by replacing E in (5.15) by q. Of course, in this way the final result depends 
on an ambiguous parameter in contrast to the energy loss of a heavy quark. For 
values of q between 5 and 20 GeV, T = 250 MeV, and a s = 0.2 the energy loss of 
a quark jet amounts to about 0.2 - 0.3 GeV/fm, agreeing well with the estimate of 
Bjorken (1982). The energy loss of a gluon jet is obtained from the one of a quark jet 
by multiplying with a color factor 9/4, as can be seen by comparing the scattering 
matrix elements (Bjorken, 1982). 

5.6 Radiative energy loss 

In the case of ultrarelativistic charged particles the radiation loss dominates over the 
collisional one. In the case of quark we then have to consider the bremsstrahlung 
diagrams in fig. 5. 9. A rigorous treatment, on the same level as for the collisional 
energy loss in the last section, of these contributions to the energy loss has not been 
succeeded so far. Gunion and Bertsch (1982) estimated the matrix elements, belong- 
ing to the diagrams in fig. 5. 9, by factorizing the process into an elastic scattering and 
a subsequent gluon emission. The energy loss per collision is not longer caused by 
the energy transfer in the elastic scattering but by the energy of the emitted gluon. 
Therefore the infrared divergence from the elastic scattering contribution is not re- 
duced but is still quadratically within naive perturbation theory as in the case of the 
damping rate of a hard quark. For a rough estimate we cut off this divergence, as 
usual, by the Debye mass. Then the order in the coupling constant is reduced by one 
order in a s , as it is also the case for the damping rate. Consequently, the radiative 
energy loss is not of order a 3 s} as naively expected from fig. 5. 9, but proportional to 
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Figure 5.9: Gluon bremsstrahlung diagrams contributing to the radiative energy loss 
of a quark in the QGP. 



a s , i.e. of the same order as the collisional energy loss. 

On the other hand, in the case of the radiative energy loss another medium effect 
besides Debye screening plays an important role. The emission of photons or gluons 
from an ultrarelativistic charged particle in a dense medium is suppressed by the finite 
formation time: A photon or gluon cannot be emitted if there is not enough time for 
the photon or gluon to reach the mass shell before the next scattering of the charged 
particle takes place. This effect predicted by Landau and Pomeranshuk (Landau and 
Pomeranshuk, 1953; Feinberg, 1966) has been verified recently experimentally in the 
transition of ultrarelativistic ions through thick targets (Klein, 1994). Since this effect 
is caused by multiple scattering events, the radiative energy loss cannot be derived 
solely from the already complicated diagrams of fig. 5. 9. First rough estimates by 
Gyulassy and Pliimer (1990) indicated that the radiative energy loss of an energetic 
parton in the QGP should be smaller than the collisional one. Later on more accurate 
investigations (Gyulassy et al., 1992; Gyulassy and Wang, 1994), however, led to a 
value for the radiative energy loss of a light quark of about 1-2 GeV/fm, suggesting 
the dominance of the radiative energy loss over the collisional one, despite of the 
presence of the Landau-Pomeranshuk effect. The radiative energy loss turned out to 
be proportional to the square of the Debye mass according to these investigations, 
which can be understood in the following way: A larger screening mass corresponds 
to a larger mean free path of the quarks in the QGP, leading to a suppression of 
the Landau-Pomeranshuk effect and thus to an enhancement of the gluon emission 
probability. Of course, if the mean free path becomes too large, the energy loss will 
decrease as bremsstrahlung processes become rare. 

In conclusion, the energy loss of a parton in a QGP is likely to be of the same 
order in a QGP and hadronic matter. Consequently, jet quenching presumably cannot 
be used gnature for the phase of the fireball. A possible exception could be 

realized for a system close to the phase transition, where lattice calculations show a 
strong reduction of the Debye mass (Gao, 1990), leading perhaps to a small energy 
loss. Then variations of the jet quenching with the temperature of the fireball or 
the center of mass energy in the heavy ion collision could hint at a possible phase 
transition (Gyulassy et al., 1992). Nevertheless, jets should be taken into account as 
a valuable external probe for the fireball in heavy ion experiments. 
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Chapter 6 

Transport rates and viscosity 



In the last chapter we have studied the effect of a relativistic plasma on an energetic 
charged particle. Now we will consider an important property of the plasma itself, 
namely its reaction to a deviation from the equilibrium momentum distribution. The 
momentum relaxation is of great interest, as it is closely related to thermalization 
times and the viscosity of the plasma. 

6.1 Motivation and definition 

The momentum relaxation follows from the collision term of the Boltzmann equation 
(Reif, 1965), which consists of a momentum integration over the differential cross sec- 
tion of the scattering process under consideration and a phase space factor containing 
the non-equilibrium distributions fi = /(r, p;,t). In the case of a two-body collision 
(1,2 <-> 3,4) this phase space factor reads (see e.g. Cassing and Mosel, 1990) 

F = h f 2 (1 ± / 3 ) (1 ± f A ) - (1 ± /i) (1 ± f 2 ) f 3 h , (6-1) 

where the upper sign applies to bosons and the lower one to fermions, respectively. 
This phase space factor vanishes not only in the equilibrium case, i.e. fi = nB,F(Pi) } 
but also if pi = p3 (or pi = P4) holds. For in this case also P2 = P4 (or P2 = P3) are 
valid due to momentum conservation, leading to f 3 = fi and / 4 = f 2 (or f 3 = f 2 and 
i*4 = fi)- Consequently, colinear elastic scattering events with scattering angle 6 = 0° 
in the center of mass system and anti-colinear with 6 = 180° do not contribute to the 
momentum relaxation. (An illuminating interpretation is given by the fact that an 
anisotropic momentum distribution cannot be equilibrated by colinear or anti-colinear 
collisions.) Since, however, those scattering angles lead to infrared singularities in 
the cross sections (6 = 0° corresponds to t = (P\ — P3) 2 = and 6 = 180° to 
u = (P1 — P4) 2 = 0), the phase space factor F softens these infrared divergences. In the 
case of a non-relativistic plasma with Coulomb interaction the cross section a can be 
replaced by a transport cross section cr trans = J da (1 — cos 6) (Lifshitz and Pitaevskii, 
1981), where the transport factor 1 — cos 6 is generated by an expansion of the collision 
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term for small scattering angles (6 = 0°), describing the suppression of large distant 
collisions. In this way the quadratic infrared singularity of the Rutherford cross 
section reduces to a logarithmic one (Coulomb logarithm) (Lifshitz and Pitaevskii, 
1981). In order to cut off also the divergence due to anti-colinear collisions, we replace 
the factor 1 — cos 6 by (sin#) 2 /2, which agree for small 6. 

In the following we will assume that momentum relaxation in a relativistic plasma 
can also be described by replacing cross sections by transport cross sections (Danie- 
lewicz and Gyulassy, 1985). This implies a modification of the interaction rate Y to 
a transport rate T trans = J dT (sin#) 2 /2. 

In the limit of small momentum transfers q = p — p', coming back to the notations 
for the momentum variables used before (see fig.4.1), the transport factor can be 
written as 2q 2 (1 - co 2 /q 2 )/s with s = (P + Kf (Thoma, 1994a). The factor of 
q 2 } appearing here, reduces the quadratic divergence in the interaction rate using 
bare propagators to a logarithmic one. Thus the transport rate turns out to be 
finite, applying the Braaten-Pisarski method analogously to the energy loss. (Owing 
to the t-u-symmetry of the relevant matrix elements and of the transport factor, 
the divergence for large momentum transfers, i.e. anti-colinear scattering events, is 
removed automatically.) The explicit computation of the transport rate can be done 
again by means of the Braaten-Yuan prescription. 

The transport rate provides an estimate for the thermalization time of the QGP in 
relativistic heavy ion collisions, T t herm = trans, because the momentum relaxation 
describes the approach to a thermal momentum distribution. This definition of the 
thermalization time has the great advantage over the one by means of the interaction 
rate, Ttherm = 1 /r, in section 4.1 that the transport rate is finite even in the absence 
of a static magnetic screening. Of course, this estimate of the thermalization time is 
justified only for small deviations from the equilibrium. On the other hand, it does 
not depend on assumptions about the begin and end of the pre-equilibrium phase, as 
it is the case using computer simulations of relativistic heavy ion collisions (Geiger 
and Miiller, 1992; Wang and Gyulassy, 1991). Furthermore, the transport rate is 
closely related to the shear viscosity of the plasma, which describes the momentum 
transport in the plasma (Reif, 1965), as it will be discussed in section 6.3. 

6.2 Transport rates 

Here we want to compute transport rates of quarks and gluons with thermal ener- 
gies, which are based on the parton scattering processes to lowest order perturbation 
theory, shown in fig. 6.1 (Combridge, Kripfganz, and Ranft, 1977; Cutler and Sivers, 
1978), where we will include screening effects again by applying an effective gluon 
propagator. (In the transport rates only the bare gluon propagator in the t- and 
u-channel diagrams leads to a logarithmic infrared singularity, while the bare quark 
propagator produces no divergence even for a vanishing quark mass.) According to 
the Braaten-Yuan prescription we decompose the transport rate into a soft (q < q*) 
and a hard (q > q*) part. 
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The soft part follows again from the quark or gluon self energy containing an ef- 
fective gluon propagator. We only have to insert the transport factor (sin#) 2 /2 under 
the integral determining the self energy. The transport factor can be approximated 
by the small momentum transfer limit 2q 2 (1 — to 2 /q 2 )/s because of q* <C t. Inserting 
this factor under the integral in (4.13), extended to QCD, we obtain 
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where d = Cf = 4/3 holds for quarks and d = Ca = 3 for gluons. The integration 
over q can be performed analytically, the one over u has to be done numerically. For 
simplicity sake we consider first only (/(/-scattering in a pure gluon gas (iVj = 0). 
Then we find (Thoma, 1994a) 
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If we are interested only in a logarithmic accuracy (Thoma, 1991a; Thoma, 1993a), 
we are able to present the final result immediately, since the hard contribution ac- 
cording to the Braaten-Yuan prescription has to assume the form 
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yielding within logarithmic accuracy 
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However, it is desirable to know the coefficient of l/a s under the logarithm, as 
a s is not really small compared to 1. For this purpose we have to compute the hard 
part, i.e. the constant Ah ar( i } explicitly. Again we start from the definition of the 
interaction rate via the matrix element containing the processes displayed in fig.6.1. 
Then e.g. the gluon transport rate in a pure gluon gas is given by 
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(6.6) 



where the matrix elements contains only the scattering diagrams (d) of fig.6.1. How- 
ever, the momentum integrations are much more involved now than in the case of the 
energy loss of a heavy quark, where we utilized M 3> T for simplifying the integrals. 
An analytic estimation can be obtained, using the approximation T trans ~ /Ocr trans , 
where p denotes the parton density of the QGP. In the case of a pure gluon gas 
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Figure 6.1: (/(/-scattering processes to lowest order. 



p = 16/ d 3 k/ (2tt) 3 riB(k) holds. In order to match the soft part onto the hard one 
(^-cancellation), the integrand of p has to be multiplied by the Bose enhancement 
factor 1 + riB(k): 

/ffil. r ( fln\ 2iii 

™lGn B {k)[l + n B (k)] Jdt(^) — . (6.7) 

Here the usual Mandelstam variables s = (P + K) 2 } t = (P — P 1 ) 2 = —2pp' (1 — cos #), 
and u = —s — t have been introduced and the transport factor has been expressed by 
these variables (sin#) 2 /2 = 2tu/s 2 (Thoma, 1994a). 

The formula (6.7) can be derived (Thoma, 1994a) setting riB(p') = in the exact 
expression (6.6), which is approximately justified because of (p 1 ) ~ 3T, choosing 
riB(k') = riB(k) } which holds approximately as long as the momentum transfer q = 
|k' — k| is not too large, and adopting the definition of the differential cross section 
(Bjorken and Drell, 1964). 

The differential cross section for (/(/-scattering has been given by Combridge, 
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Kripfganz, and Ranft (1977): 

9g 4 f us st tu 




Combining (6.7) and (6.8) and using the limits —s-\-q* 2 and —q* 2 lor the t-integration 
results in 
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(6.9) 



where we made use of s 3> q* 2 . In the case of thermal gluons we replace s by its 
thermal average (s) = 2(p)(k) = 14.59T 2 . After adding the soft contribution (6.3), 
we end up with the q* independent transport rate 

25 

T g , t rans(N f = 0) ~ 5.2 a 2 T In— . (6.10) 

a s 

The computation of the gluon and quark transport rates in a QGP with Nj = 2 
is a little bit more elaborate, leading to (Thoma, 1994a) 

19 

T g ,trans(Nf=2) ~ 6.6 a 2 s T In — — , 

a s 

21 

T q , t rans(N f = 2) ~ 2.5 a 2 s T In — — . (6.11) 

a s 

These results for the transport rates of thermal partons, (6.10) and (6.11), have 
the unpleasant property of being negative if extrapolated to realistic values of the 
coupling constant (a s = 0.2 - 0.5). For energetic particles, however, with s 3> T 2 
the transport rates remain positive (see (6.9)). This behavior, demonstrating the 
shortcoming of the logarithmic approximation, is caused by isolating the leading 
order contribution by means of the Braaten-Yuan prescription, analogously to the 
unphysical results for the energy loss of a heavy quark (section 5.5) for small momenta. 
It can be avoided only by going beyond the leading order a 2 s} in which case we have 
to consider effective vertices in addition. 

Similar results for the transport rates in the QGP have been found by Heisel- 
berg (1994a), becoming unphysical if a s > 0.3. Heiselberg derived the transport 
rates numerically from a variational solution of the Boltzmann equation taking into 
account an effective gluon propagator in the matrix element of the collision term. 
However, as he did not use the Braaten-Yuan prescription but integrated over the 
entire momentum range of the exchanged gluon, it is questionable if his final result is 
gauge independent and free of higher order contributions. In any case, his elaborate 
consideration of the transport theory by means of a variational method justifies the 
introduction of a transport factor, as his results shows the same dependence on the 
coupling constant and the temperature and the same infrared behavior (Heiselberg 
and Pethick, 1993) as in (6.10) and (6.11), with which it also agrees quantitatively 
quite well. 
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Figure 6.2: Velocity gradient in a plasma leading to the shear viscosity. 

In conclusion, the Braaten-Yuan prescription allows for a calculation beyond the 
leading logarithm, if one is capable to evaluate the hard part. Whereas we succeeded 
to derive the hard contribution in the case of the energy loss of a heavy quark exactly, 
we had to adopt some approximations for the hard part of the transport rate. Anyway, 
the final result, coinciding with the one by Heiselberg (1994a), indicates that the 
leading order calculations for the transport rates are inconsistent, if extrapolated to 
realistic values of the strong coupling constant. Therefore, one has to be careful using 
results, such as thermalization times, obtained only within logarithmic accuracy. 

6.3 Viscosity 

Let us consider a plasma that deviates from the equilibrium state by a velocity gradi- 
ent as in fig. 6. 2. The shear viscosity coefficient t] indicates how fast the plasma equi- 
librates due to momentum transport between the different velocity domains. This 
behavior can be investigated by means of the transport theory. From elementary 
kinetic considerations or using the relaxation time approximation for the Boltzmann 
equation 

m = 3 Y.f><- (p^ A « ( 6 - 12 ) 

i 

has been found in the non-relativistic case (Reif, 1965), where pi denotes the particle 
density, (pi) the mean momentum, and A; the mean free path of the particles of kind 
i. Solving the Boltzmann equation for the deviation from the equilibrium, described 
above, by a variational calculations or expansion in moments (Reif, 1965), leads to a 
modification of the formula (6.12) by introducing the transport factor (sin#) 2 in the 
cross section entering the mean free path and by replacing the factor of 1/3 by 0.327. 
This formula gives a result, for instance, in the case of a gas of hard spheres which 
deviates only by 1.6% from the exact one, whereas (6.12) underestimates the exact 
result by a factor of 1.5 (Reif, 1965). 

In a relativistic plasma the factor of 1/3 in (6.12) has to be replaced by 4/15 (de 
Groot, van Leeuwen, and van Weert, 1980; Danielewicz and Gyulassy, 1985). The 
mean free path A; follows from the transport rate A; = 1/ (2Ti jtrans ) } where the factor 
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Figure 6.3: Shear viscosity coefficient of the QGP as function of the strong coupling 
constant according to (6.13) (solid line) and Baym et al. (1990) (dashed line). The 
Navier-Stokes equation is valid below the dotted line. 

of 1/2 is caused by the choice of the transport factor (sin#) 2 /2 instead of (sin#) 2 in 



Inserting the quark and gluon transport rates (6.11) in A; in (6.12), we find for 
the shear viscosity coefficient of the QGP (Thoma, 1994a) 



The contribution from the quarks to the shear viscosity is about three times that of 
the gluons, since gluons interact stronger, having therefore a shorter mean free path. 
In fig. 6. 3 rj/T 3 is depicted as a function of the coupling constant (solid line). The 
dashed line corresponds to a calculation by Baym et al. (1990), which is based on a 
variational solution of the Boltzmann equation. However, in this calculation the vis- 
cosity has been determined only up to logarithmic accuracy; i.e. the coefficient of 1 /a s 
under the logarithm has been set equal to 1, resulting in r\ = 1.16T 3 /[a 2 s ln(l/a s )]. 
Consequently, both the calculations in fig. 6. 3 deviate significantly from each other for 
a s > 0.1, where the complete result (6.13) breaks down. 

Similar results within logarithmic accuracy of the same form as the one of Baym 
et al. (1990) have been reported in the literature with coefficients 0.28 (Hosoya and 
Kajantie, 1985), of 0.57 (Danielewicz and Gyulassy, 1985), and 1.02 (Thoma, 1991a) 
instead of 1.16. The first two results have been obtained using the Debye mass as an 
infrared regulator, whereas the latter result relies on the Braaten-Pisarski method. 
The dependence on the coupling constant, i] ~ 1/a 2 , up to a logarithmic factor can 
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be traced back to the use of the transport rate (see also Mrowczyhski (1990)). Using 
instead the interaction rate for the mean free path, rj ~ l/a s , is found. In other words, 
the simple formula (6.12) without transport factor leads to a result that deviates by 
a order of a s from the one improved by the transport factor. 

An alternative access to the viscosity (Hosoya, Sakagami, and Takao, 1984; Furu- 
sawa et al., 1984; Horsely and Schoenmaker, 1987; Ilyin et al., 1989) besides the Boltz- 
mann equation is provided by the Kubo formulas, in which the transport coefficients 
can be derived from equilibrium correlation functions of the energy-momentum tensor 
in accordance with the fluctuation-dissipation theorem (Zubarev, 1974). This method 
first led to a result with a wrong dependence on the coupling constant, rj > 2.6T 3 j a s 
(Ilyin, Panferov, and Sinukov, 1989), that later on has been corrected (Ilyin et al, 
1992). This original discrepancy between both the approaches was the reason for the 
investigation of the viscosity within the Braaten-Pisarski method (Thoma, 1991a), 
which allows for a consistent treatment of the infrared behavior which is responsible 
for the dependence on the coupling constant. 

In fig. 6. 3 the limit of the validity of the Navier-Stokes equation, providing a hydro- 
dynamic description of the expansion phase of the QGP in ultrarelativistic heavy ion 
collisions, is shown in addition. As opposed to the Euler equation, the Navier-Stokes 
equations considers dissipative effects such as the viscosity. According to Danielewicz 
and Gyulassy (1985) the Navier-Stokes equations, not to mention the Euler equation, 
can be used only if the viscosity is smaller than rj < er/4, where e is the energy density 
of the QGP and r the proper time. Since r = (t /Tq)T 3 holds, neglecting dissipative 
effects (Bjorken, 1983), we get the upper limit, shown in fig. 6. 3, for optimistically 
chosen values of To = 1 fm/c and T = 1 GeV. The perturbatively derived results are 
above this limit even in the case of logarithmic accuracy. This means that hydrody- 
namical calculations starting from an ideal fluid without dissipation (Euler equation) 
are questionable. (The consideration of dissipation (Navier-Stokes equation) in the 
case ultrarelativistic heavy ion collisions seems to be difficult (Strottman, 1994).) 

Also simple models taking account of non-perturbative effects just above the 
critical temperature, T = 1 - 2 T c , led to an estimate for the viscosity, which shows a 
non-negligible entropy production, and which should be considered in hydrodynamic 
models (Cleymans et al., 1994). Lattice calculations of the viscosity, based on the 
Kubo formulas, are affected by large uncertainties so far (Karsch and Wyld, 1987). 
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Chapter 7 

Damping and production of hard 
photons 

Electromagnetic probes (photons, dileptons) have been proposed as promising sig- 
natures of the QGP formation (Ruuskanen, 1992). Photons and dileptons produced 
in relativistic heavy ion collisions leave the fireball without further interactions, thus 
providing a direct probe for the strongly interacting system. As an example we will 
consider the thermal production of hard photons with energies E 3> T and momen- 
tum p from a QGP. To lowest order perturbation theory it is caused by the Compton 
scattering like process of fig. 7. 1(a) and the quark pair annihilation of fig. 7. 1(b). Ne- 
glecting the quark masses the production rate E dR/d 3 p suffers from a logarithmic in- 
frared divergence. However, this singularity is screened by medium effects, leading to 
an effective temperature dependent quark mass of the order of gT (see (2.22)), which 
cuts off the divergence much more effectively than the bare mass. This medium effect 
can be taken into account consistently using the effective quark propagator (2.29). 

The damping rate *y(p) of a hard photon is defined as the imaginary part of the 
dispersion relation of a real, i.e. transverse, photon, *y(p) = —Imuj(p) } which follows 
from the transverse part of the photon self energy, containing a quark loop. To lowest 
order perturbation theory the mechanism of photon damping is given by the inverse 
inverse processes of fig. 7.1, i.e. Compton scattering and pair creation with a photon 
in the initial state. These processes are related to the definition of the damping rate 
via the photon self energy, shown in fig. 7. 2, by cutting rules, similarly as in the case 
of the damping rates of hard leptons and partons, discussed in section 4.1. (Cutting 
the self energy in fig. 7. 2 leads to diagrams describing the production as well as the 
damping of a photon in the QGP (Weldon, 1983).) 

Photon production and damping are related to each other by the principle of 
detailed balance (Weldon, 1983): EdR/d 3 p = [4/(2vr) 3 ] exp(-£/T) 7 (p), where the 
factor of 4/(27r) 3 is a matter of definition. As a matter of fact, the calculation of 
the damping rate is easier than the one of the production rate, as discussed below. 
Hence we will present only the calculation of the damping rate here, following Thoma 
(1994b). 
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Figure 7.1: Photon production in the QGP to lowest order. 



7.1 Photon damping rate 



Since the photon damping rate to lowest order naive perturbation theory is loga- 
rithmically infrared divergent for a vanishing quark mass, a hnite result complete 
to leading order in the coupling constants can be obtained again by means of the 
Braaten-Yuan prescription. Then the soft contribution follows from using the effec- 
tive quark propagator in the photon self energy and the hard one from the matrix 
elements of the inverse processes of fig. 7.1. Using the Braaten-Pisarski method for the 
soft part, the photon self energy of fig. 7. 3 will be considered, where the blob denotes 
the effective quark propagator. In the case of hard photons one effective quark prop- 
agator is sufficient. Also no effective quark-photon vertex is needed. The imaginary 
part of the photon self energy according to fig. 7. 3 comes from the imaginary part of 
the effective quark propagator, corresponding to the exchange of a collective quark 
mode. 

The dispersion relation of a hard real photon results from 



where the transverse part of the photon self energy is given in (2.15). From (7.1) we 
get for the damping rate of a hard photon in the case of no overdamping, 7 <C p, 



scale in the exchanged quark momentum K , follows from using the photon self energy 
shown in fig. 7. 3 in (7.2), given by 



pl-p 2 - n T (po,p) = , 






(7.3) 
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Figure 7.2: Lowest order contribution of the photon self energy in naive perturbation 
theory to the photon production and damping in the QGP. 



where the factor of 10/3 comes from adding up both the diagrams in fig. 7. 3, from 
summing over the colors of the internal quarks, and considering the electric charges of 
up and down quarks, assumed to be thermalized in the QGP. The quark propagators 
S(Q = P — K) and S*(K) denote the bare and the effective one, respectively, for 
which we will adopt the helicity representation of (2.29). 

Evaluating the trace over the 7 matrices, we get for the transverse part of the 
polarization tensor 



10 , f d 4 K 
n T = — i e 2 



3 J 2tt 



1 fl-V 1 + V\ 1 fl + V 1-V 



D+(K) \d+(Q) d_(Q)J D-(K) \d+(Q) d_(Q) 

where the functions D±(K) are given in (2.30), d±(Q) = —qo ± q, and V = (p • k)(p • 
q)/(p 2 qk). 

Using the imaginary time formalism, the integration over k in (7.4) is replaced 
by a discrete sum, which can be evaluated most easily by adopting the spectral repre- 
sentation for the fermion propagators (Braaten, Pisarski, and Yuan, 1990; Kapusta, 
Lichard, and Seibert, 1991), analogously to the one of the photon propagator, (4.6), 



1 r ~dTe koT l~ duj p±(w,k)[l-n F (ui)] 



D±(K 



1 1-0 foo 

dre qoT / du'r±(u>',q) [1 -n F (w')] T , (7.5) 



d±(Q) 

where we only need the discontinuous part of the spectral function, p dts (uj,k) = 
j3±(uj 7 k)6(k 2 — io 2 ) with f3± = —Im{l/ D±)/ir, as the damping is due to the imaginary 
part of the effective quark propagator, and r±(a/, q) = &[w' ± q). Using furthermore 
p > T and < k 2 - to 2 < k* 2 < T 2 , we find 

Imn T (po,p) = --£- I dk [ duj[(k-uj)/3 + (uj } k) + (k + uj)/3_(uj } k)] } (7.6) 



The integral in (7.6) also appears in the soft part of the photon production rate 
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Figure 7.3: Photon self energy with an effective quark propagator determining the 
soft part of the photon damping. 

(Kapusta, Lichard, and Seibert, 1991). It can be evaluated numerically, leading to 



5e 2 , , k* 2 
Isojt = — rn 



4 In 

2i7rp q 2m 



q 

5vr aa s T 2 , 0.3183&* 2 

= T^ ln ^r^' (7 - 7) 

where the effective quark mass is given by m 2 = g 2 T 2 /6. 

The hard contribution, for which the momentum K of the exchanged quark is 
larger than /.•' . can be calculated most conveniently from the matrix elements accord- 
ing to the inverse processes of fig. 7.1, which are Compton scattering 

(2tt) 4 6\P + K -P'- K') 12 {\M\ 2 ) comp (7.8) 
and pair creation 

(2tt) 4 6\P + K - P> - A") 16 (\M\ 2 ) pair . (7.9) 

Here momenta with a prime belong to outgoing particles. The factors in front of 
the amplitudes come from summing over the possible states of the incoming thermal 
parton. The amplitudes, averaged over initial states and summed over final ones, 
using Mandelstam variables, s = (P + K) , t = (P — P') 2 , u = —s — t, are given by 

(\M\ 2 )comp = -— f-V (- + - 

21 \s u 

(\M\% mr = \e 2 g 2 + , (7.10) 
9 \t uj 

where the amplitude of each process is the sum of the one involving an up quark and 
a down quark. (These contributions have to be added incoherently, as the electric 
charges of these quarks are different.) The dependence on the Mandelstam variables 
can be taken over from the corresponding QED processes (Halzen and Martin, 1984), 
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modified by color factors from averaging and summing over the initial and final color 
states, respectively. 

In order to evaluate the integrations over the final states in (7.8) and (7.9), we 
assume p 1 ^> T and k' 3> T, which is justified since p > p' + k' 3> T and the phase 
space for p'~T or &'~T is unfavorable (Kapusta, Lichard, and Seibert, 1991). Thus we 
may assume 1 + rag (p') ~ l — npi^p') ~ 1 — npik 1 ) ~ 1 corresponding to the Boltzmann 
approximation for the distribution functions, e.g. 1 + riB(p') = exp(p'/T) riB(p') — 1, 
also used by Kapusta, Lichard, and Seibert (1991) in the case of photon production. 
This approximation simplifies the expressions (7.8) and (7.9) considerably, because 
we can now evaluate the integrals over p' and k' transforming to the center of mass 
system, using the Lorentz invariant phase space factor (Halzen and Martin, 1984) 

dL = (2tt) 4 6\P + K -P'- K') , d Y , d l k ' = — . (7.11) 
V ; V ; (27r) 3 2p' (27r) 3 2k> 8tts y ' 

In the calculation of the photon production rate one has to integrate over the ini- 
tial momenta. These integrations are rather involved (Kapusta, Lichard, and Seibert, 
1991; Staadt, Greiner, and Rafelski, 1986) due to the distribution functions of the 
incoming particles, even assuming the Boltzmann approximation. However, in the 
case of photon damping, where one integrates over the final momenta (see (7.8) and 
(7.9)), the distribution functions of the outgoing particles vanish within the Boltz- 
mann approximation. (There is no Pauli blocking or Bose enhancement for classical 
particles described by Boltzmann distributions.) Hence, these integrations can be 
performed easily using (7.11). 

Let us first consider Compton scattering: 

dL(\M\ 2 ) Comp = ^e 2 g 2 (m-^ + £| , (7.12) 

where k* 2 <C s cuts off the logarithmic infrared divergence of the it-channel. (In order 
to match the soft part onto the hard part consistently, we have chosen a covariant 
separation scale in both the contributions, which facilitates the calculation of the 
hard part.) Using s = 2pk(l — p ■ k) we obtain, after integrating over d 3 k } 



Comp _ 57T aa s T 2 



^fhard 



27 p 



\ ^pT 1 C'(2) 



(7.13) 



where 7 = 0.57722 is Euler's constant and ((z) is Riemann's zeta function with 
C'(2)/((2) = -0.56996 (Braaten and Thoma, 1991a). 
The pair creation contribution gives analogously 

IOtt aa s T 2 \ 4pT ('(2) 
In — r - 1 - 7 + 



pair 
lhard 



27 p 



k* 2 1 C(2) 



(7.14) 



Adding up the (7.13) and (7.14), we get the total hard contribution 



5tt aa s T 2 n 0.9706pT 
lhard = ~. In — 2 . 7.15 

9 p k* 
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Adding up the soft contribution (7.7) and the hard one (7.15), the separation 
scale k* cancels, as expected, yielding the final result (Thoma, 1994b) 



5tt aa s T 2 0.2317p 
7 = - In — . 7.16 

9 p a s l 

Again the extrapolation to realistic values of a s fails if the photon momentum becomes 
too small, p < 4.4a s T, as a consequence of isolating the leading order term by the 
Braaten-Yuan prescription. 

The mean free path, A = I/7, of a 1 GeV photon in a QGP at T = 200 MeV, 
extrapolating the result (7.16), obtained in the weak coupling limit j < 1, to a 
realistic a s = 0.3, is A = 480 fm. This value is much larger than the dimensions of a 
fireball in ultrarelativistic heavy ion collisions, confirming that photons may be used 
as a direct probe of the fireball (Shuryak, 1978). (Of course, higher order processes 
such as bremsstrahlung may contribute to the rate similarly for realistic values of the 
strong coupling constant.) 



7.2 Photon production rate 

As mentioned above, the production rate of a hard photon in the QGP is related to 
the damping rate by the principle of detailed balance. In this way we obtain from 
(7.16) 

„dR 5 2 _ E/T 0.2317£ 

E ir = T^ aasTe ln — ^~ ' ( 7 - 17 ) 

to lowest order in the coupling constants, denoting the photon energy by E instead 
of p now. The same result has been derived by Kapusta, Lichard, and Seibert (1991) 
and by Baier et al. (1992) directly. Their calculation, however, is somewhat more 
tedious, since they had to integrate over the Boltzmann distributions of the initial 
states, where they could not make use of (7.11). Kapusta, Lichard, and Seibert 
(1991) also justified the Boltzmann approximation by using the exact distributions, 
integrating numerically. 

Kapusta, Lichard, and Seibert (1991) also estimated the thermal photon pro- 
duction rate from hadronic matter, concluding that the QGP shines as brightly as 
hadronic matter at the same temperature. Recent investigations, including additional 
mesons in the hadron gas (Xiong, Shuryak, and Brown, 1992; Haglin, 1994), indicate 
that the hadronic systems emits at least twice as many photons as the QGP at the 
same temperature. 

These rates can be used in hydrodynamical models of the evolution of the fireball, 
assuming a QGP as well as a hadronic phase, for deriving photon spectra, measurable 
in heavy ion experiments. In this way a good agreement with the data from SPS has 
been found already (Arbex et al., 1994). 

However, these estimates have not yet included pre-equilibrium effects. Although 
thermalization, i.e. isotropic and exponential parton distribution functions, seems to 
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be achieved rapidly, already after about 0.2 - 0.3 fm/c, the chemical equilibration 
takes place much later, if at all (Biro et al., 1993; Geiger, 1993). On one side, an 
early thermalization corresponds to an high initial temperature (600 MeV at RHIC, 
800 at LHC (Biro et al., 1993)), increasing the photon rate, on the other side, an 
incomplete chemical equilibration, in particular for the quark component, leads to 
a suppression of the photon rate. Phenomenological investigations of the influence 
of these pre-equilibrium effects show that the emission of high-energy photons are 
dominated by the pre-equilibrium phase (Kampfer and Pavlenko, 1993; Strickland, 
1994). Similar results hold for dileptons (Kampfer and Pavlenko, 1992; Kawrakow 
and Ranft, 1992). A non-equilibrium resummation taking into account medium effects 
analogously to the Braaten-Pisarski method, allowing for a consistent calculation of 
the pre-equilibrium photon production, has not been developed so far. 

Furthermore, effects of a finite quark chemical potential should be considered for 
AGS-, SPS-, and maybe also RHIC-energies, as the nuclear transparency is expected 
to be incomplete for these energies. For example, the RQMD code (Keitz et al., 1991) 
indicates that the quark chemical potential even in a QGP formed at RHIC might be 
of the order // = 1 - 2T close to the mid-rapidity region (Dumitru et al., 1992). Since 
lattice calculations do not allow for considering a finite chemical potential up to now, 
the only way of doing this is given by perturbative methods. 

The Braaten-Pisarski method can be extended to finite chemical potential in a 
straightforward way (Vija and Thoma, 1994). The only modification of the effective 
Green's functions is due to the effective masses, which now depends besides the tem- 
perature also on the chemical potential (Braaten and Pisarski, 1992a). The effective 
gluon mass is given by (Toimela, 1985; Toimela, 1986) 



where ///is the quark chemical potential of quarks with flavor /. (In the following 



In addition, the Fermi distributions, for example used in the hard parts, deter- 
mined by integrations over matrix elements, also contain now the chemical potential: 
np(E p ± fi) = l/[exp((£ , p ± fi)/T) + 1], where the minus sign indicates quarks and 
the plus sign anti-quarks. 

In this way the hard photon production rate has also been derived for a QGP with 
a finite quark chemical potential /t (Traxler, Vija, and Thoma, 1994). The soft part 
of the photon production (or damping) rate at finite quark chemical potential follows 
immediately from the // = case by substituting for m q in (7.7) the expression (7.19). 
In order to match the hard part onto the soft one (^-cancellation), we have to give up 






The effective quark mass reads (Kajantie and Ruuskanen, 1983) 




(7.19) 
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the Boltzmann approximation in the hard part, as it produces a //-independent pre- 
factor in front of the logarithm, as e.g. in (7.15). Using exact distribution functions 
for all particles, e.g. in (7.8) and (7.9), the integrations in the hard part have to be 
done numerically. The final result can be fitted by (Traxler, Vija, and Thoma, 1994) 

w dR 5 T2 (, . I" 2 \ -e/t , 0-2317£ 

E ^- p = 18^ aasT { 1 + ^) e ln ^^' (7 - 20) 

which is independent of the separation scale as expected. Surprisingly also the 
chemical potential drops out under the logarithm within an error of about 3%. This 
cancellation is due to adding the soft and the hard part, as the soft part contains // 
under the logarithm according to (7.7). So the chemical potential shows up only in 
front of the logarithm, coming from the pre-factor m 2 q in the soft part, (7.7). 

If we put into (7.20) typical numbers for the quark chemical potential expected 
at RHIC, fi = 1 - 2T, the modification due to the chemical potential is negligible 
(fig. 7.4), as the photon production rate is dominated by the exponential decrease 
with the temperature and the effect of the chemical potential is suppressed compared 
to the temperature by a factor of 1/tt 2 according to (7.20). However, if we regard 
the energy density, e = e(T, //), instead of the temperature as given, since it can be 
derived more directly from the center of mass energy y/s of the heavy ion collision, 
e ~ \n.y/s (Satz, 1992), than the temperature (Dumitru et al., 1993a), a significant 
reduction of the photon production rate results (fig. 7. 5). This suppression is caused 
by the fact that an increase of the quark chemical potential reduces the temperature 
for a given energy density, as can be seen for example by adopting a MIT-bag model 
equation of state (Dumitru, 1993b), 

f37ir2 U7raA T 4 + 3 (l - ^) TV + A (l " — V + B (7.21) 



V 30 3 J V 7T J " 27T 2 V 7T 

with a bag constant B = (200 MeV) 4 . In conclusion, the temperature is the domi- 
nating factor compared with the chemical potential. 

The production rate of soft photons has also been investigated by means of the 
Braaten-Pisarski method (Baier, Peigne, and Schiff, 1994). Here effective vertices as 
well as propagators have to be used. However, in this case an infrared divergence, 
arising from a mass singularity in the effective photon vertex, has been encountered. 
This problem has not been solved so far. 

Also the production of soft dileptons has been calculated, using the effective per- 
turbation theory (Braaten, Pisarski, and Yuan, 1990; Wong, 1992). The calculation 
is similar to the one of the soft photon production, but there is no infrared singularity, 
as the virtual photons decaying into dileptons are not on mass shell. An interesting 
structure of the soft dilepton spectrum has been found that originates from the dis- 
persion relations of the collective quark modes. Unfortunately, this structure cannot 
be used presumably for detecting the QGP, since it is covered by a large background. 
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Figure 7.4: Photon production rate of the QGP according to (7.20) as a function of 
the photon energy at a given temperature T = 200 MeV, a s = 0.4, and various values 
of the quark chemical potential //. 
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Figure 7.5: Photon production rate of the QGP according to (7.20) as a function of 
the photon energy at a given energy density e = 5 GeV/fm 3 , a s = 0.4, and various 
values of the quark chemical potential //. 
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Chapter 8 
Conclusions 



The development and application of field theoretic methods at finite temperature 
(and chemical potential) has been pushed recently by present (AGS, SPS) and future 
(RHIC, LHC) heavy ion experiments at relativistic energies (y/s = 5 - 6300 GeV/A). 
In particular, these theoretical investigations are of crucial importance for the search 
for the QGP (Miiller, 1993). 

Perturbative QCD can be extended to thermal systems by considering imaginary 
times and discrete energies (Matsubara frequencies) (Kapusta, 1989). Compared to 
lattice calculations this method has the advantage of being capable to treat dynamical 
properties of the QGP, such as the most signatures proposed. For the evaluation of 
loop integrals within the imaginary time formalism it is most convenient to adopt a 
mixed time-momentum representation of the propagators (Saclay method) (Pisarski, 
1988). Important examples for the application of these perturbative methods are the 
self energies of plasma particles (Klimov, 1982; Weldon, 1982a, 1982b), from which 
dispersion relations, damping and production rates follow. 

The use of perturbative methods for gauge theories, however, is complicated by 
gauge symmetry and massless gauge bosons: Naive perturbation theory, i.e. the 
exclusive use of bare propagators and vertices, may lead to gauge dependent and 
infrared divergent results for physical quantities. A famous example is given by the 
gluon damping rate at rest (plasmon puzzle) (Lopez, Parikh, and Siemens, 1985). 

Braaten and Pisarski (1990a) first recognized that these problems of the pertur- 
bation theory are based on the fact that at finite temperature diagrams of higher 
order may contribute to lower order in the coupling constant. In other words, results 
gained by using bare Green's functions can be incomplete in the order of the coupling 
constant. In order to avoid this problem, one has to resum a certain class of infinite 
many diagrams, the HTL, which are characterized by hard loop momenta of the order 
T, corresponding to the high temperature limit. This resummation can be performed 
by means of the Dyson-Schwinger equation and leads to effective propagators and ver- 
tices, showing a complex momentum dependence. These effective Green's functions 
can be used in a perturbative expansion analogously to bare ones. However, this is 
necessary only if all external momenta of the Green's function under consideration 
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can be soft, i.e. of the order of gT . 

The resummation of certain diagrams for avoiding infrared singularities is known 
for a long time in the case of non-relativistic problems and the scalar field theory. 
The achievements of Braaten and Pisarski consist in extracting the relevant diagrams 
(HTL) for gauge theories, which is much more complicated due to the momentum 
dependence of the HTL and Ward identities than for the scalar field theory, and in 
proving that the effective perturbation theory produces gauge independent results for 
observables (Braaten and Pisarski, 1990c). 

Applying the Braaten-Pisarski method it is possible to obtain consistent results, 
which are gauge independent and complete to leading order in the coupling constant, 
and which show an improved infrared behavior. (Many quantities that are infrared 
divergent within naive perturbation theory are finite now.) The physical meaning 
of using effective Green's functions is the inclusion of medium effects of the QGP, 
such as Debye screening by the color charges of the thermal partons. In this way 
Braaten and Pisarski (1990b) succeeded in deriving the gluon damping rate at zero 
momentum consistently, thus solving the plasmon puzzle. Thereby they found that 
the damping mechanism is completely different from the one suggested in naive per- 
turbation theory, based on g gg processes. Using effective propagators and vertices 
for computing the imaginary part of the gluon self energy, from which the damping 
rate follows, it can be seen by cutting the diagrams of the self energy that gg <->■ gg 
and gg <->■ ggg processes are responsible for the damping. The imaginary part of the 
gluon polarisation tensor arises from the imaginary part of the effective Green's func- 
tions, which in turn comes from another medium effect, caused by the resummation, 
namely Landau damping. 

In conclusion, the Braaten-Pisarski method means an enormous progress for the 
qualitative understanding of the QGP as well as explicit calculations compared to 
naive perturbation theory. The application of the effective perturbation theory to 
realistic situations of the QGP, however, is hampered by two difficulties. 

First, the Braaten-Pisarski method is based on the assumption of the weak cou- 
pling limit, j < 1, allowing for a distinction between hard and soft momenta, whereas 
realistic values are given by g = 1.5 - 2.5. Thus the extrapolation of results obtained 
by using effective Green's functions to realistic values of g is questionable. On the 
other hand, the assumption, ^ < 1, might be too restrictive and should be replaced 
by a s <C 7r (Braaten, 1995). This could be checked by considering corrections of the 
next to leading order in g. So far only the plasma frequency, i.e. the thermal photon 
mass, has been calculated beyond the leading order. In the case of a pure gluon gas it 
is given by = m 2 (1 — 0.31g) with m 2 = g 2 T 2 /3 (Schulz, 1994), implying that the 
correction for g = 1 is about 30%. Furthermore, quantities that are logarithmically 
infrared divergent using bare propagators but finite applying the Braaten-Pisarski 
method, such as the energy loss, the transport rate, or the photon production rate, 
become unphysical, i.e. negative, above g ~ 1, but only if the energies of the particles 
under consideration are not much larger than the temperature. This observation in- 
dicates that results achieved by means of the effective perturbation theory within the 
weak coupling limit can be extrapolated to realistic values of the coupling constant, 
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at least lor quantities ol energetic particles with E 3> T. For thermal particles, how- 
ever, g < 1 has to be required (Thoma, 1994d). This restriction, lor example, renders 
the determination ol thermalization times in ultrarelativistic heavy ion collisions by 
means ol the transport rate impossible, indicating the lailure ol the leading logarithm 
approximation in this case. Further hints lor the validity ol the extrapolation ol 
the results to realistic values ol the coupling constant can be found by comparing 
with non-perturbative results, such as lattice calculations and the mean field ap- 
proximation. For example, the equation ol state for the QGP obtained Irom lattice 
QCD coincides with the one derived Irom perturbative QCD above about twice to 
four times the critical temperature (Petersson, 1991). Also the effective gluon mass, 
found within the Hartree approximation (Thoma, 1992b, 1993b), differs Irom the 
perturbative result to lowest order by less than 30% at g = 1.5. Summarizing, the 
Braaten-Pisarski method not only considers important medium effects in the QGP 
qualitatively but also allows for explicit calculations in many cases, leading to results 
which are quantitatively correct within a lactor ol, let me say, about 2. 

The second problem, in contrast to the one discussed above, is a conceptual one, 
also occurring in the weak coupling limit, namely the absence ol a static magnetic 
screening in the effective gluon propagator. Therefore quantities which are sensi- 
tive to the momentum scale g 2 T are still inlrared divergent even at leading order, 
using effective Green's lunctions. For example the damping rates ol hard partons 
are quadratically inlrared divergent in naive perturbation theory and logarithmi- 
cally in the effective one. Assuming an inlrared cut-off ol the order g 2 T } e.g. a 
non-perturbative magnetic screening mass (Pisarski, 1989b) or a finite width ol the 
partons in the QGP (Lebedev and Smilga, 1990), the Braaten-Pisarski method al- 
lows for a result ol these rates within logarithmic accuracy. Up to now only one 
observable, the color conductivity (Selikhov and Gyulassy, 1993; Heiselberg, 1994b), 
has been found, which suffers Irom this problem to leading order. However, Irom 
a certain order on probably all quantities become sensitive to the order g 2 T } where 
the Braaten-Pisarski method is insufficient. To what extent these higher order cor- 
rections are important alter regularization by a method beyond the Braaten-Pisarski 
method is an open question. Such a method might be based on a sell-consistent 
resummation, such as the Hartree approximation (Thoma, 1992b, 1993b). However, 
the Hartree approximation leads to inconsistent gauge dependent results and can be 
regarded only as a starting point for more elaborate approximations, ol which only 
preliminary investigations exist (Biro, 1989; Kalashnikov, 1992; Mishra et al., 1993; 
Simonov, 1993). 

Further non-perturbative effects ol the QGP, which cannot be described by the 
Braaten-Pisarski method, play a role at and just above the phase transition, such 
as the survival ol mesonic and baryonic correlations in the deconfined phase (Boyd, 
1994; Koch et al., 1992). Furthermore, there are predictions ol non-perturbative 
field configurations, such as instantons (Gross, Pisarski, and Yaffe, 1981), monopoles 
(Oleszczuck and Polonyi, 1992; Biro and Miiller, 1993), or glue balls (Rischke et al., 
1992) even at arbitrary high temperatures. Within the Braaten-Pisarski method we 
are able to describe important deviations Irom an ideal gas, e.g. by medium effects 
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like Debye screening, at least. Indeed, there are indications that the effective gluon 
mass may explain deviation from an ideal gas, observed on the lattice (Biro, 1992; 
Goloviznin and Satz, 1993; Peshier et al., 1994). 

Finally, the Braaten-Pisarski method can be applied only to equilibrated systems. 
A complete equilibration of the QGP in relativistic heavy ion collisions, however, 
may never happen (Biro et al., 1993). Non-equilibrium situations can also be treated 
diagramatically by means of the Keldysh formalism (Lifshitz and Pitaevskii, 1981). 
A resummation of certain diagrams analogously to the Braaten-Pisarski technique 
might be imaginable (Altherr, 1994). 

Thermal quantities, which diverge only logarithmically in naive perturbation the- 
ory, can be calculated consistently beyond the logarithmic approximation within the 
Braaten-Pisarski method in the weak coupling limit. This can be achieved by using 
the Braaten-Yuan prescription (Braaten and Yuan, 1991), where the soft part for mo- 
mentum transfers below a separation scale q* is evaluated using effective propagators, 
whereas for the hard part above q* bare Green's functions are sufficient. Restricting 
the arbitrary separation scale by gT <C q* <C T, it cancels after adding up both the 
contributions, leading to a result which is complete to leading order in the coupling 
constant. 

A basic quantity, which has to be treated by means of the Braaten-Pisarski 
method and which represents the starting point for important observables of the 
QGP, is the damping (or interaction) rate of a quark or gluon. These rates can be 
derived either from the self energy of the particles under consideration or from the 
matrix elements of the corresponding processes. In the case of hard particles with 
momenta at least of the order of the temperature, where an effective gluon propaga- 
tor suffices, a logarithmically infrared divergence due to the missing static magnetic 
screening has been found (Thoma and Gyulassy, 1991). Moreover, these rates are 
proportional to g 2 instead of g 4 } as expected in naive perturbation theory. This is 
caused by the appearance of the square of the effective gluon mass, m 2 ~ g 2 T 2 } in 
the denominator of the effective gluon propagator and the fact that soft momenta 
contribute significantly in the integrals defining the rates. However, it is doubtful 
whether the damping rates have a physical relevance by themselves (Lebedev and 
Smilga, 1990). 

Damping rates of partons at zero momenta are infrared finite. They require the 
use of effective propagators and vertices at the same time, rendering their calculation 
very tedious (Braaten and Pisarski, 1990b, 1992b). 

The collisional energy loss of an energetic parton in the QGP (Bjorken, 1982; 
Thoma and Gyulassy, 1991), caused by elastic scattering via the exchange of a gluon, 
follows from the hard interaction rate by inserting the energy transfer per collision, 
divided by the velocity of the parton, under the integral defining the rate. In this 
way a field theoretic definition of the energy loss has been achieved, extending the 
Bethe-Bloch formula to the energy loss in a relativistic plasma (Braaten and Thoma, 
1991a). 

The energy loss of high-energy partons in a QGP is of interest as it represents 
an important input parameter for the evolution of jets in ultrarelativistic heavy ion 
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collisions. In particular, preliminary estimates indicated that the energy loss in the 
QGP might be significantly smaller than the one in a hadron gas, providing therefore 
a signature for the QGP ("jet unquenching" ) (Gyulassy and Pliimer, 1990). This 
speculation, however, has been revoked after including the radiative energy loss due 
to bremsstrahlung, which was estimated to be of the same order as the energy loss in 
hadronic matter (Gyulassy et al., 1992). 

Using the Braaten-Yuan prescription, the collisional energy loss of a heavy quark 
in the QGP has been calculated consistently (Braaten and Thoma, 1991b); i.e. a 
finite, gauge independent result has been derived, complete to leading order in the 
coupling constant and free of ambiguous parameters, such as a maximum or minimum 
momentum transfer or an arbitrary separation scale. Owing to the additional factor 
of the energy transfer per collision compared to the interaction rate, the energy loss 
is not only finite, using an effective gluon propagator in the soft part, but also of 
the order of g 4 in contrast to the interaction rate. The independent evaluation of 
the soft and the hard part according to the Braaten-Yuan prescription allows for a 
result beyond the logarithmic approximation. For small energies below a few GeV the 
extrapolation to realistic values of g breaks down. There are, however, no problems 
for energies above about 10 GeV. The collisional energy loss of light partons requires 
the introduction of an empirical maximum energy transfer per collision, in order to 
describe the jet evolution, yielding typical values of dE / dx = 0.3 GeV/fm (Thoma, 
1991b). 

Another interesting quantity of the QGP, which can be derived by means of 
the Braaten-Yuan prescription in a similar way as the energy loss, is the transport 
rate, from which thermalization times and the shear viscosity follow (Thoma, 1991a, 
1994a). It can be obtained from the interaction rate by introducing a transport 
factor, which suppresses colinear and anti-colinear scattering events in the center of 
mass system, describing momentum relaxation. The transport factor improves the 
infrared behavior of the interaction rate in the same way as the energy transfer per 
collision in the case of the energy loss. Therefore a finite expression of the order 
g 4 results. Unfortunately, it is not possible to obtain a meaningful estimate of the 
thermalization time from this, as the transport rates become negative if extrapolated 
to realistic values of the coupling constant for thermal energies. 

The shear viscosity coefficient, which is inversely proportional to the transport 
rate, has been calculated in this way beyond the logarithmic approximation (Thoma, 
1994a). The result indicates that dissipative effects in the expansion phase of the 
QGP are important (Danielewicz and Gyulassy, 1985). 

Further interesting quantities, calculable consistently by means of the Braaten- 
Yuan prescription, are the production and damping rates of hard photons in the 
QGP (Kapusta, Lichard, and Seibert, 1991; Baier et al., 1992; Thoma, 1994b). The 
damping rate determines the mean free path of the photons in the QGP, while the 
production rate describes the thermal photon emission of the QGP, which might 
serve as a signature for the QGP. The soft part now requires the use of an effective 
quark propagator instead of a gluon propagator, as for the energy loss or transport 
rate, because the damping (production) of hard photons is caused to lowest order by 



73 



Compton scattering and pair creation (annihilation). As a result the mean free path 
of a photon in the QGP is confirmed to be much larger than the dimensions of the 
fireball (Thoma, 1994b), showing that hard photons provide a direct probe for the 
fireball. On the other hand, a hadron gas seems to emit about as much photons as 
a QGP at the same temperature (Kapusta, Lichard, and Seibert, 1991), suggesting 
that the photon production can be used only as a thermometer of the fireball but not 
for distinguishing between the deconfined and the hadronic phase. 

The introduction of a finite quark chemical potential, describing a surplus of 
quarks over anti-quarks in the QGP, originating from an incomplete transparency of 
the nuclear collision, into the Braaten-Pisarski method is straightforward (Vija and 
Thoma, 1994). The effective Green's functions are only changed by new effective 
gluon and quark masses, now containing also the chemical potential besides the tem- 
perature. The computation of the hard parts, however, is more complicated due to 
the integration over Fermi distributions with finite chemical potential. In this way 
the energy loss (Vija and Thoma, 1994) and the photon production rate (Traxler, 
Vija, and Thoma, 1994) have been extended to a finite quark chemical potential. 
Only minor modifications due to the chemical potential, even for // = 2T, have been 
observed, demonstrating that the temperature plays the dominant role. If we keep, 
however, the energy density instead of the temperature fixed, the temperature de- 
creases with increasing chemical potential, leading to a significant reduction of the 
results, in particular in the case of the photon production (Dumitru et al., 1993b). 
Which quantity should be considered as given, the temperature or the energy density, 
has to be decided by the experiment. 

A large number of interesting quantities of a QGP, but also of a QED or elec- 
troweak plasma, which are relevant for astrophysical problems (see e.g. Altherr and 
Kraemmer (1992) and Braaten (1992)), remain to be calculated by using the Braaten- 
Pisarski method. In particular, a systematic consideration of medium effects in the 
elementary parton cross sections is desirable. In this way, for example, the production 
and diffusion of strange or charm quarks could be treated consistently. 

Furthermore, it would be of great interest to extend the Braaten-Pisarski re- 
summation to the pre-equilibrium phase, which might dominate many signatures in 
the fireball, for example by starting from a diagramatic non-equilibrium technique 
(Keldysh formalism). Finally, the development of methods beyond the Braaten- 
Pisarski method, suitable for describing non-perturbative phenomena (magnetic gluon 
mass, hadronic correlators, instantons, etc.), maybe within a self-consistent resum- 
mation, is still missing. 
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Appendix A 
Notations 



Here the notations used in this report are summarized. In addition, some Feynman 
rules are presented. 

Natural units, i.e. % = c = 1, are adopted throughout the paper. Furthermore, 
Boltzmann's constant has been chosen as k = 1. 

The metric tensor is given by the Minkowski metric g^ with the diagonal elements 
goo = 1 and g lt = g 22 = #33 = -1. 

Greek indices z/, ... of covariant vectors denote the space-time components, Latin 
ones z, j, ... only space components. 

Four momenta are indicated by capital letters, the magnitudes of three momenta 
by small letters, e.g. K 2 = k^ — k 2 with k = |k|. This notation is convenient at finite 
temperature, as magnitudes of three momenta appear often, e.g. in the distribution 
functions riB,F(k). 

The bare propagator of a scalar particle with mass m is depicted by a dashed line 
and given by 



The bare propagator of a fermion with mass M , represented by a solid line with 
an arrow, reads 



The bare propagator of a photon is indicated by a wavy line, the one of a gluon 
by a spiral line. The propagators of gauge bosons depend on the gauge choice and 
are given when needed (see e.g. (2.25)). 

The bare vertex of the ^-theory is given by i 4! g 2 and the photon-electron vertex 
by i e^f^. The vertices of QCD have not been used explicitly. They can be found e.g. 
in Itzykson and Zuber (1980). 

In the case of the ^-theory symmetry factors have to be considered, e.g. a factor 
of 1/2 for the tadpole in (2.7). 

Otherwise standard Feynman rules are used (see e.g. Mandl and Shaw (1984)). 
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Appendix B 

Calculation of the photon self 
energy 



Here the longitudinal part of the photon self energy in the HTL-limit is calculated 
explicitly, because it is one of the most important quantities of gauge theories at finite 
temperature. The computation exemplifies also the perturbative methods at finite 
temperature, developed in chapter 2. We follow the ideas of Braaten and Pisarski 
(1990a) mainly. 

The starting point is given by the diagram in fig. 2. 2. Using standard Feynman 
rules, first at zero temperature, we get 

M(P) = -*e 2 l^tr [ 7 , S(K - P) 7 „ S(K)] . (B.l) 



At finite temperature, after evaluating the trace over the 7 matrices and using (2.11), 
we obtain 

00 r d 3 k ~ 
IV(P) = -4 e 2 T £ j^(K,Q l , + Q,K l ,-g^K-Q)A(K)A(Q), (B.2) 

n= — 00 V ^/ 

where Q = P — K and k and q assume discrete values, e.g. k = (2n + l)iirT . 

Restricting ourselves to the longitudinal component, il o = 11^, and introducing 
the Saclay propagators, (2.12), we find 

n L (P) = -4e 2 /^r^Mo + k-q) [ P dTdT'e k ° T e^'A(T,k)A(T',q) , (B.3) 

where we neglected the fermion mass (T 3> m e ). 

In order to evaluate the sum over k 0} we absorb the term proportional to k q in 
derivatives according to r and r' and integrate by parts, 

Il L (P) = -4 e 2 [^L T J2l P dr dr> e^-') 
J {2ttY ^ Jo 



k 

dA(r,k) dA{r',q) . ~ , 7\T/ / \ 
j \ , + k • q A(r, k) A(r', q) 



;b.4) 
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After summing over k by means of 



T J2 e ko{T ~ T > = S(t-t') (B.5) 

Tl— — CO 

and integrating over r ! afterwards, we get using (2.13) 

X {n F (k) [1 - n F (q)] e {k - q)T + [1 - n F (k)] n F (q) e^ k+q)T ) . (B.6) 

Here we assumed already k • q/ '(kq) ~ —1, which follows from the high temperature 
limit, as discussed below. This assumption leads to a cancellation of half of the terms 
in the sum of both the propagator products in (B.4), since the derivative according 
to r changes the sign between the terms in (2.13). 
The integration over r yields 

(i , , e P(±k±q) _ 1 

dr e { P° ±k± ^ T = - , (B.7) 

p ± k ± q 

where we utilized p = 27rmT, as the particle belonging to p is a boson. Now the 
summation over k is concluded and we may continue the remaining p analytically; 
i.e., from now on p is regarded as real, continuous energy of the photon. By means 
of the useful relation 

nF(k) e~? k , (B.8) 



1 — n F (k) 



(B.6) reduces to 



H L (P) = 2 e 2 / -^y 3 [n F (k) - n F (q)] (—± L— ) . (B.9) 

J (Ziry \p + k-q p - k + q J 

This integral cannot be evaluated analytically. In the high temperature limit, T ^> p 0} 
p } however, the integral can be simplified. For this purpose we decompose the integral 
by introducing a separation momentum k* } restricted by p <^ k* <^ T . Owing to 
dimensional reasons, the integral from to k* is proportional to k* . The other 
possible scale in the integral, namely the temperature, drops out of the calculation 
because of n F (k) ~ n F (q) ~ 1/2 if k } q <^ T. Then the contribution from k < k* can 
be neglected, since the final result in the high temperature limit, (2.16), is of order 
T 2 , as we will see. For the integral over the hard momenta, k > k* } the approximation 
q ~ k — pi] with j] = p ■ k can be adopted, from which we find 

1. k • q/ (kq) ~ —1, as used already in (B.6), 

2. p ± (k — q) ~ p ± prj } and 

3. n F (k) — n F (q) ~ prj dn F (k) / dk. 

It should be noted that the expression (B.9) contains no temperature independent 
term. This is only the case, because we have used there already the approximation 1. 
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On the other hand, the temperature independent term has to be proportional to P 2 
after renormalization and can be neglected, too. Next we substitute the approxima- 
tions 2. and 3. into (B.9) and integrate over k from zero to infinity. The error made 
in this way is negligible, as it is of order k* 2 . Consequently, the high temperature 
limit, T 3> p, po } is equivalent to integrating over hard loop momenta (HTL), k ~ T. 
After integrating over rj and k, using 



2 dn F (k) tt 2 T 2 



dk k 2 „ v ' = , (B.10) 

dk 6 y ' 



we end up with, (2.16), i.e. 



n I (P ) = -f!f(i-fii n 2!i±£), (B.n; 

3 V 2p po-pj 
from which we see that the integral in (B.9) is of order T 2 , indeed. 
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